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Abstract. A local moment approach is developed for single-particle excitations of 
a symmetric Anderson impurity model (AIM) with a soft-gap hybridization vanishing 
at the Fermi level: Ai cx , with r > 0. Local moments are introduced explicitly 
from the outset, and a two-self-energy description is employed in which single-particle 
excitations are coupled dynamically to low-energy transverse spin fluctuations. The 
resultant theory is applicable on all energy scales, and captures both the spin- 
fluctuation regime of strong coupling (large-[/), as well as the weak coupling regime 
where it is perturbatively exact for those r-domains in which perturbation theory 
in U is non-singular. While the primary emphasis is on single-particle dynamics, the 
quantum phase transition between strong coupling (SC) and local moment (LM) phases 
can also be addressed directly; for the spin-fluctuation regime in particular a number 
of asymptotically exact results are thereby obtained, notably for the behaviour of the 
critical Uc{r) separating SC/LM states and the Kondo scale u)ra{r) characteristic of the 
SC phase. Results for both single-particle spectra and SC/LM phase boundaries are 
found to agree well with recent numerical renormalization group (NRG) studies; and a 
number of further testable predictions are made. Single-particle spectra are examined 
systematically for both SC and LM states; in particular, for all Q < r < i, spectra 
characteristic of the SC state are predicted to exhibit an r-dependent universal scaling 
form as the SC/LM phase boundary is approached and the Kondo scale vanishes. 
Results for the 'normal' r = AIM are moreover recovered smoothly from the limit 
r — > 0, where the resultant description of single-particle dynamics includes recovery of 
Doniach-Sunjic tails in the wings of the Kondo resonance, as well as characteristic low- 
energy Fermi liquid behaviour and the exponential diminution with U of the Kondo 
scale itself. The normal AIM is found to represent a particular case of more generic 
behaviour characteristic of the r > SC phase which, in agreement with conclusions 
drawn from recent NRG work, may be viewed as a non-trivial but natural generalization 
of Fermi liquid physics. 
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1. Introduction 

The Anderson impurity model (AIM) [1] is the archetype for describing dilute, correlated 
magnetic impurities in metals. Reviewed comprehensively in [2], its essential strong 
coupling behaviour is that of the Kondo effect: the spin-| impurity is quenched by 
coupling to low-energy excitations of the non-interacting metallic host. But while 
thermodynamic properties of the model are well understood, the same cannot be 
said for a theoretical description of dynamics, in particular those of single-particle 
excitations. Here a wide variety of theories have been developed, including the non- 
crossing approximation (NCA) [3-6], expansions [7-9] and slave boson approaches 
[10-12]. Their undoubted successes notwithstanding, however, each has significant 
limitations. They are designed to capture the iV — > oo limit (as opposed to the spin-i 
case, N = 2), and the extreme asymmetric limit of C/ = oo with U the local (impurity) 
interaction: extension to finite-C/ is not straightforward. The NCA describes well high- 
energy single-particle excitations embodied in the Hubbard satellite bands, but fails to 
capture the low-energy Fermi liquid behaviour of the Kondo resonance. Slave boson 
approaches by contrast are much less satisfactory on high-energy scales, their virtue 
being to handle low-energies and Fermi liquid behaviour. Even here however they are 
limited: in common with 1/A^ expansions, they cannot for example recover the Doniach- 
Sunjic tails [13] in the Kondo resonance that arc observed in numerical renormalization 
group [14] and quantum Monte Carlo [15] studies of the AIM, and known to be important 
experimentally [16]. 

The normal AIM has of course one 'simplifying' feature: the host is metaUic by 
presumption, whence the host-impurity coupling embodied in the hybridization function 
A(a;) is essentially frequency-independent and controlled by its value at the Fermi level, 
u; = 0; in consequence, and excepting the atomic limit where the impurity/host trivially 
decouple, the impurity spin is quenched and the system is a conventional Fermi liquid 
for all f/ > (see e.g. [2]). The brief remarks above are nonetheless not confined to this 
problem, but symptomatic in general of the widely accepted need for new theoretical 
approaches to strongly correlated electrons. And the challenge is naturally more acute 
when A(ci;) acquires a frequency dependence that can lead to qualitatively new physics, 
and in particular the possibility of nontrivial zero-temperature phase transitions. One 
example of such arises in the lattice-based Hubbard model, which within the framework 
of Dynamical Mean-Field Theory (or the infinite-dimensional limit) maps onto an 
effective AIM [17] with a hybridization that is a functional of the impurity Green 
function, and is hence both cj-dependent and to be determined self-consistently. The 
quantum phase transition here is the celebrated and still controversial Mott transition, 
occurring at a critical 11^ in the paramagnetic phase of the half-filled (particle- hole 
symmetric) model: between a gapless Fermi liquid metal, and a gapped local moment 
insulator characterized by an /cBln2 residual entropy (for reviews, see e.g. [18-20]). 

A second example, considered in this paper, is the soft-gap AIM appropriate to 
a semi-metallic host, or one that itself may be viewed as being on the verge of a 
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simple band-crossing metal-insulator transition: in which the (imaginary part of the) 
hybridization function exhibits a soft-gap at the Fermi level, Ai(cj) oc \ujY with r > 0, 
in contrast to the normal 'metallic' AIM, r = 0, for which Ai(c(j = 0) is constant. This 
problem, for which a wide range of possible experimental candidates arise (see eg [21]), 
was first studied by Withoff and Fradkin [22] in the context of the corresponding soft- 
gap Kendo model, using both 'poor man's' scahng and a large-A^ mean-field theory. 
Much study of the soft-gap Kondo and Anderson models has since ensued, in particular 
via scahng [22,23], large- expansions [22,24,25], the numerical renormalization group 
(NRG) [21,23,26-28] and perturbation theory in U [29]. It is known thereby that two 
distinct ground states exist, between which in general a quantum phase transition occurs 
at a finite critical Uc{r): a doubly degenerate local moment (LM) state in which the 
impurity spin remains unquenched; and a strong coupling (SC) state in which the 
impurity spin is locally quenched, and a Kondo cfi'cct is manifest. 

The underlying physics is known to be particularly rich for the particle-hole 
symmetric model, to which NRG studies in particular have devoted considerable 
attention, including both thermodynamic properties [21,26-28] and impurity single- 
particle spectra [28]. It is the symmetric spin-| soft-gap AIM that we consider here, 
by developing a microscopic 'local moment approach' that has recently been applied 
successfully to the normal (r = 0) AIM [30]. Our primary focus is thus an analytical 
treatment of single-particle dynamics — on all energy scales, and for any interaction 
strength U — as embodied in the impurity Green function, G{uj), and hence spectrum 
D{uj) oc ImG(u;); although an integral clement of the approach also permits statics, in 
the form of the SC/LM transition and associated phase boundaries, to be addressed 
directly. 

The resultant theory, which seeks in particular to capture the spin-fiuctuation 
regime of strong coupling (large-C/) but is also perturbatively exact in weak coupling for 
those r-domains in which straight perturbation theory in U is known to be applicable 
[29], appears to be rather successful. Its predictions for single-particle spectra and 
SC/LM phase boundaries agree well, both qualitatively and quantitatively, with extant 
NRG calculations [21,28] that provide essentially exact numerical results for the problem. 
In the spin-fluctuation regime of strong coupling a number of asymptotically exact 
results are obtained, in particular for the Kondo scale characteristic of the SC phase 
and the low-r behaviour of the critical Uc{r). Results for the normal AIM [30] are 
moreover obtained smoothly from the limit r — > 0, where the resultant description 
of single-particle dynamics includes recovery of the Doniach-Sunjic tails, as well as 
characteristic Iow-cl^ Fermi liquid behaviour and the exponentiality of the Kondo scale. 
We find in fact that the normal AIM constitutes in many ways just a particular case 
of more generic behaviour characteristic of the r > SC state, which we argue may 
be regarded as a 'generalized Fermi liquid', in agreement with the conclusions drawn 
by Gonzalez- Buxton and Ingersent [21] from detailed NRG studies. One manifestation 
of this is the prediction that, for all r e [0, |), single-particle spectra characteristic of 
the SC state should acquire a universal scaling form as the SC/LM phase boundary is 
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approached, thus generahzing to finite-r behaviour that is famihar in the context of the 
normal AIM [14,15]; this and related predictions of the present theory will be tested 
against NRG calculations in a subsequent publication [31]. 

The paper is organized as follows. A brief introduction to the soft-gap AIM is 
given in section 2, where two facets are highlighted. First, the non-interacting limit 
U — Its behaviour, in contrast to that of the normal AIM, is nontrivial: both SC 
and LM states arise, for r < 1 and r > 1 respectively, and with distinct signatures in 
the underlying spectral functions [29]. Second, we emphasise the generalized pinning 
condition estabhshed by us previously [29], whereby A{ijj) — \ijjYD{ijj) is pinned at the 
Fermi level cu = for any r and U where a SC state obtains; and which represents a 
generalization of the corresponding condition familiar for the r = AIM where it is 
normally viewed as a consequence of the Friedel sum rule (see e.g. [2]). Imposition of 
this spectral pinning as a self-consistency condition plays a central role in the current 
work. 

In section 3 we introduce the 'two-self-energy' description that underlies the present 
theory. Such an approach is physically natural if one aims to describe the doubly 
degenerate LM state; and, we would argue, is at least desirable if one seeks to construct 
a non-perturbative theory that can simultaneously handle the possibility of both LM 
and SC states, and hence the transition between them. As for the normal AIM [30], 
our approach to the interaction self-energies starts from the simplest non-trivial mean- 
field approximation in which the notion of an impurity local moment, determined self- 
consistently, is introduced explicitly from the outset: unrestricted Hartrcc-Fock (UHF), 
as considered by Anderson in his original paper [1] on the r = AIM. The deficiencies 
of this static mean-field approximation per se are of course severe; but it is in large part 
a physical understanding of them, considered in section 4.2, that enables a subsequent 
many-body approach to be developed successfully. Moreover, and in contrast to the 
normal AIM, even UHF is nontrivial for the soft-gap problem: as shown in section 4.1 
it gives rise for example to a phase diagram that, in predicting for all finite-f/ solely LM 
states for any r > |, concurs qualitatively with the results of NRG calculations [21,28]. 

Dynamical many-body contributions to the self-energies, over and above the Fock 
term alone retained at simple mean-field level, are detailed in section 5 for both SC and 
LM states. In physical terms these embody coupling of single-particle excitations to 
low-energy transverse spin fluctuations, and capture the dynamical spin-flip scattering 
required in particular to describe the Kondo, or spin-fluctuation, regime. SC states are 
obtained via self-consistent imposition of the generalized pinning condition, the limits of 
stability of solutions to which in turn give the critical Uc{r) for the SC/LM transition. 
In this way, as for the r = AIM [30], the Kondo or spin-flip scale ou^ = u}^{r) arises 
naturally within the present approach, and is found to be non-zero throughout the entire 
SC phase, vanishing continuously only as the SC/LM phase boundary is approached, 
U — » Uc{r)—. For U > Uc{r) in the LM phase, by contrast, ci;in(r) = 0: as expected 
physically for a doubly degenerate state, where there is no energy cost for a local spin 
flip. 
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Resultant phase boundaries between SC and LM states are considered in section 6; 
including (section 6.3) their predicted one-parameter scaling behaviour in the regimes 
U/D ^ 1 and ^ 1 (with D the bandwidth of the host spectrum or hybridization 
function), and detailed comparison with extant NRG results [21,28] (section 6.4). The 
evolution and critical behaviour of the SC Kondo scale u!in{r) is considered in section 6.1. 
Particular emphasis is given here to small r -C 1 where salient results can be extracted 
analytically; and which is of evident importance in connecting to the normal r — AIM, 
for which the exact exponential asymptotics of the Kondo scale are correctly recovered. 
As r — > 0, the critical Uc{r) for the SC/LM boundary is found to be oc 1/r, and the 
condition for the Schricffcr- Wolff transformation [32] mapping the soft-gap AIM to the 
corresponding Kondo model is thus satisfied. The critical exchange coupling Jc{r) for 
the soft-gap Kondo model as r — can thus be obtained from the present approach 
(section 6.2); and is found to be given precisely by the scaling result obtained originally 
by Withoff and Fradkin [22], which we argue is asymptotically exact. 

Single-particle impurity spectra, and their evolution with interaction strength U 
from strong to weak coupling, are considered explicitly in sections 7 and 8. The 'bare' 
D{u;) are discussed in section 7, on all energy scales and for both SC and LM states. 
Many-body broadening of the high-energy Hubbard satellites, whose existence is well 
known for the normal AIM (see eg [6,30]), is argued to arise also in the soft-gap problem 
and shown to be correctly recovered by the present approach; as too are the characteristic 
u spectral signatures of the SC and LM phases found in NRG calculations [28], viz 
D{u}) ~ and ~ |(x;|'' respectively. The behaviour of the spectra in weak coupling, 

U —>■ 0, is considered in section 7.1. For r < | where both the ^7 — > and the ^7 = 
ground states are SC, and for r > 1 where the ground state is found to be LM for 
all [/ > 0, the theory is shown to be perturbatively exact to (and including) second- 
order in U about the noninteracting limit. For | < r < 1 by contrast, the U > 
ground state is found to be LM but the non-interacting ground state is SC [29]; and 
the natural breakdown of conventional perturbation theory in U is clearly evident in 
the non-analyticity (in U) of the conventional 'single' self-energy as f/ — 0. Explicit 
comparison to single-particle spectra obtained from NRG calculations [28] is made in 
section 7.2, and excellent agreement is found. 

Finally, we consider in section 8 the spectral functions A{u) = |ci;|''£'(ci;) in the SC 
phase — 'modified' to remove the \uj\^^' divergence at low-o; that is symptomatic of the 
SC state, and which is entirely unrenormalized by interaction effects [29]. The A{u;) are 
found to exhibit familiar characteristics: they are pinned at the Fermi level, a; = 0, and 
contain a generalized Kondo resonance whose width is proportional to the Kondo scale 
u!ra{r) and thus narrows progressively as U is increased towards the SC/LM phase where 
uj^{r) vanishes. This is just as found for the normal AIM (see eg [2,14,15]) where the 
critical (f//Ao)r=o = oo (reflecting the fact that the SC/LM transition here coincides 
trivially with the atomic limit). Moreover, as for the r = AIM, we find generally 
for any given r < | that as U ^ Uc{r) — , the A{u!) becomes a universal function of 
a;/a;m(r); with a scaled Kondo resonance that exhibits characteristic r-dependent low- 
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frequency behaviour, as well as Doniach-Sunjic tails in the 'wings' of the spectrum for 
^/^m(r) > 1. 

2. Background 

We begin with some necessary background, particularly in relation to the non-interacting 
limit (section 2.1), and the generalized spectral pinning condition characteristic of a SC 
state with U >0 [29] (section 2.2). 

With the Fermi level taken as the energy origin, the Hamiltonian for the spin-| 
Anderson model is given in standard notation by 

H ^^host ~l~ -^impurity ~l~ -^^hybridization 

= ^ ^kflka + ^{^i + \Uni-a)nia + ^ Vik{c\^CUc + C^kaCia) (2.1) 

with efe the host dispersion, Vik the hybridization and the impurity level; for the 
symmetric Anderson model considered here, = —U/2 with U the on-site interaction. 
We focus on the zero-temperature single-particle impurity Green function, defined 

by 

G{t) = -i(T{Q.(t)cL}) = G+(t) + G-{t) (2.2) 

and separated for later purposes into retarded and advanced components; since H is 
invariant under a —a, G is naturally independent of spin, a. We add that while 
the primary physical content of G{u) is that of single particle dynamics, analysis of it 
will also enable identification of the phase boundaries between SC and LM states, as 
detailed in sections 5 and 6. 

2.1. Non-interacting limit 

For U — 0, the impurity Green function g{uj) is given by 

g{u) = [u + \r]sgn{uj) - A{uj)y^ : rj ^ 0+ (2.3) 
and is determined by the hybridization function 

A(a;) = T ^^"^^ . , = Ar(u;) - isgn(u;) Ai(u;) (2.4a) 

^ a; - efe ir;sgn(a;) 

with 

Ai{u;)^7rY,\VikMu;-ek). (2.4b) 

k 

We consider a symmetric hybridization, A(a;) = — A(— a;), and in particular a power-law 
form 

Ai{lo) = Ao\uj\''e{D - \uj\) (2.5) 

with r > and bandwidth D {6{x) being the unit step function). A pure power-law 
hybridization, while naturally not realistic on arbitrary scales, captures the requisite low- 
LJ behaviour in the simplest way; moreover, as familiar from the usual r = Anderson 
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model (see e.g. [2]), one expects impurity properties to be controlled primarily by 
the low-cj behaviour and largely independent of detailed band structure. Note also from 
equation (2.4b) that to specify Ai(c<j), the {V^fc} and host eigenvalues {cfc} do not require 
separate specification; but that for the particular case of constant Vik = V (considered 
in section 6.2), Ai(ti;) and the host spectrum phost('-^) are simply related: 

Ai(a;) =7rl^Vhost(a;). (2.6) 

The real part of the hybridization, A^{u!), follows from a Hilbert transform, viz 

F^iuj) = F,{co,) P (^—] (2.7) 

with F — A; and hence 

Ar(cu) = sgn(a;)Ao|a;r- / dy ^ (2.8) 

TT Jo (1 - y ) 

where a principal value is henceforth understood. We shall require explicitly only the 
low-o; behaviour of A^{u!), given from equation (2.8) (for any r > 0) by 



Ar(cu) = -sgn(cu)Ao <j tan (f r) \u;\' + — + O 



I a; I 

'd 



(2.9) 



7r(r -1) D 

Notice from equation (2.8) that the wide-band limit D — > oo, as commonly employed 
for the normal Anderson model r — (see e.g. [2]), can be taken for r < 1; and for this 
case, Aji{u;) — — sgn(a;)Aotan(|r)|a;|'' V a;. 

Given the hybridization function, the non-interacting g{ou) — Keg{u}) — 
i7rsgn{uj) doluj) follows directly from equation (2.3). All relevant impurity properties are 
determined by the spectral density (Iq^uj), including 'excess' thermodynamic functions 
induced by addition of the impurity, and local properties such as the impurity 
susceptibility XuC^) — ~9I^B{d{Siz) /dh) \h=o (with h a magnetic field acting solely 
on the impurity). Details are given in [29]; here we summarize results relevant to the 
present work. 

The key feature of the non-interacting hmit is that LM states occur exclusively for 
r > 1 and SC states for r < 1; with clear signatures of the respective phases apparent 
in the single particle spectrum do{u;). For r > 1 (LM) do{u;) is given for < D by 

do{uj) = qS{uj) + d^{uj) : r > 1 (2.10a) 

and contains both a discrete state at a; = 0, with poleweight q = [1 — {dA-[i/du)^=o]~^ 
given from equation (2.9) by 

= l + (2.10b) 
7r(r — 1) 

and a continuum (or 'band') piece (io(a;) ~ ["^T^^ as a; — > 0. The pole contribution to 
do{u;) is the characteristic spectral signature of the U — LM state. It produces [29] for 
example a local susceptibility xliT) = (g72)xcurie(7') as T ^ 0, i.e. liuiT-^oTxliT) = 
(f{g^-Qf'/'&: the impurity spin remains unquenched, symptomatic of a LM state. For 
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r < 1 (SC) by contrast, there is no a; = pole contribution and dQ{uj) = (£^{00) is given 
as a; ^ by 

d»= . r,^r2,. ,1 +Q(l^r"^0 :r<l (2.11) 
ttAo [1 + tan^(|r)J 

for any r > 0, with a characteristic divergence for < r < 1. In consequence 

[29], limT^oTxu{T) = 0: the spin is quenched as occurs for the normal r — Anderson 
model, one reason why the SC state may be regarded [21,29] as a natural generalization 
of conventional Fermi liquid physics. 

2.2. SC state: spectral pinning 

For U > the impurity Green function, G{uj) = X{u!) — msgn{u!)D{u!) {— —G{—uj) by 
particle- hole symmetry), may be expressed as 

G{uj) = [u; + ir]SgTL{uj) - A(u;) - S(u;)]"^ (2.12) 

where 

E(a;) = ER(a;) - isgn(o;)Si(o;) (2.13) 

is an interaction self-energy whose real/imaginary parts are related by the Hilbert 
transform equation (2.7) with F — Ti. (Note that S(cc;) is defined to exclude the trivial 
Hartree contribution: from particle-hole symmetry the Fermi level remains fixed at 
a; = for all [/ > (and any r), whence the impurity charge rij = X](t(^«o-) = 1 V C/ 
and thus the Hartree contribution of {U/2)ni — U/2 trivially cancels = —U/2.) 

In [29] we have established conditions upon S(a') for a SC state to arise for U > 0. 
These are, very simply, that Si(a;) (and hence ER,(a;)) should decay to zero as a; — > 
more rapidly than |a;|^ (with r < 1), i.e. 

Ei(u;) a\u\^ : A > r (2.14) 

with A > r. In consequence, as follows directly from equation (2.12), the low-frequency 
behaviour of D(oi)) is that of do(oi)); hence from equation (2.11), D(oi)) \^\~'^ — 
which is indeed the spectral signature of the SC state found in finite- C/ NRG studies 
[28]. Further, using 

lim \(jj\''D(u;) = lim |a;|%(a;) (2.15) 
and defining the modified spectral function 

A{uj) = \uj\''D{uj) (2.16a) 
equation (2.11) yields the pinning condition 

ttAo [1 + tan^d r)] A{uj = 0) = 1 (2.16b) 

for all U and r where a SC state obtains. This result will prove central to our analysis 
in the following sections. It encompasses as a special case the well known result for the 
normal r — Anderson model (see e.g. [2]): that 7rAo-D(a; = 0) = 1 — the impurity 
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spectrum is pinned at the Fermi level c<j = for any U where a normal Fermi liquid 
state obtains (in that case all [/ > 0). Equation (2.16) generalizes the pinning condition 
to arbitrary r for a SC state, and reflects the fact that interactions have no influence 
in renormalizing the asymptotic behaviour of D{uj) as ^ 0, again consistent with 
the view [21,29] that the SC state constitutes a natural generalization of Fermi liquid 
behaviour. 



3. Two-self-energy description 

Equation (2.12) merely defines a single self-energy E(a;), via a Dyson equation 

G{uj)^g{uj) + g{uj)T.{uj)G{uj) (3.1) 

that does not by itself enable a calculation of the impurity Green function. And while 
at first sight it may invite a perturbative treatment in U about the non-interacting limit 
(where G(uj) = g{oj)), there are two reasons to be wary of such an approach. First, the 
general applicability of such is not obvious in the soft-gap problem: indeed, as discussed 
in [29], there is evidence to suggest that for r G [|, 1], perturbation theory in U is 
inapplicable as [/ ^ (a point to which we return again in section 7.1). Second, and 
more generally, even if a perturbative approach in U is possible for sufficiently low U 
— as it is known to be [29] for < r < | and r > 1 — straight perturbation theory 
in the interaction strength naturally cannot capture the transition between SC and LM 
states. For this, one requires an inherently non-perturbative approach that is capable 
of describing both the doubly degenerate LM state and the SC (or generalized Fermi 
liquid) state, and hence the quantum phase transition between them. 

To this end, and noting that direct calculation of a single self-energy is not 
sacrosanct, we follow recent work [30] on the normal r = Anderson model and adopt 
a two-self-energy description, with the impurity Green function expressed formally as 

G{uj) = \[G^{uj) + G^{uj)] (3.2) 

where 

G^{uj)^ uj + ir]Sg-a{uj)- i^{uj)-t^{uj) ^ (3.3) 

with interaction self-energies Ea(a;) (and a =] / [ or +/—)■ Such a description is 
self-evidently natural to describe a doubly degenerate LM state and, more generally, is 
appropriate if one seeks to construct a non-perturbative theory starting from either the 
atomic fimit (Vik — 0, where Eo-(a;) = —\aU) or an unrestricted Hartree Fock (UHF) 
mean-field approach. It is the latter strategy that we adopt and, although the deficiences 
of a simple UHF approach by itself are significant, use of it as a starting point for a 
genuine many-body treatment will be shown to yield a rather successful description of 
the soft-gap problem, as well as the normal r = Anderson model considered hitherto 
[30]. 

For the symmetric case under consideration, particle-hole symmetry implies 

Gt(a;) = -G^{-uj) (3.4) 
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(and thus D^{(jj) = Di{—u) for the corresponding spectral densities D^{ijj) — 
— 7r~^sgn(a;)ImG(^(a;)); from which, since A(a;) = — A(— a;), equation (3.3) imphes 

t^{u) = -t^{-uj). (3.5) 

In consequence, from equations (3.2) and (2.12), the impurity Green function and single 
self-energy satisfy the familiar conditions 

^(a;) = -G{-uj), S(o^) = -E(-o^). (3.6) 

Equations (3.4-6) merely express a basic symmetry, which must of course be satisfied by 
any approximate theory; and equation (3.5) in particular shows that it is sufficient to 
consider only one of the So- (a;), say S|(tc;). Once a theory for S|(cc;) has been developed, 
direct comparison of equations (3.2,3) with equation (2.12) permits, if desired, S(a;) to 
be determined, via 



1 

2 


E|(c<j) — S|(— cij) 4- 


- 2(7(a;)E|(cj)E|(— cu) 


1 - \g{uj) 









= — p z (3.7) 

where g{uj) is the non-interacting Green function, equation (2.3). 

The are obviously not calculable exactly, but diagramatic perturbation 

theory based upon a UHF mean-field state can be employed to develop suitable (and 
indeed asymptotically exact) approximations as detailed in section 5. To this end it is 
helpful to separate the full interaction self-energies as 

E^(u;) = -|[/|//|+Eo(u;) (3.8) 

where ±|[/|/x| is the purely static Fock contribution which alone survives at UHF 
level (with the local moment magnitude); and where the Ti„{uj) — to which the 
symmetry equation (3.5) also applies — contains the dynamics that, at low frequencies 
in particular, are naturally central to the problem. 

Before proceeding we note that conditions upon Eo-(a;) for a SC state to arise for 
f/ > are readily established, independently of any specific approximation. As in section 
2.2 we consider explicitly r < 1 (since a SC state must be pcrturbatively continuable 
from the non-interacting limit and for r > 1 the U = Q ground state is a LM one [29]). 
The requisite conditions upon Eo-(a;) are identical to those of section 2.2 for E(a;): the 
real and imaginary parts of Eo-(a;) must decay to zero as a; — > more rapidly than 
Ai/R ~ \ujY. From equations (3.2,3) the a; ^ behaviour of G^{uj), Gi{uj) and G{ijj) 
then coincide, and reduce to that of the non-interacting limit; equation (2.15) is thus 
satisfied, and in consequence the generalized spectral pinning condition characteristic 
of the SC state, equation (2.16b), follows. 

Finally, notice that a necessary condition for a SC state to arise is thus Eo-(cj = 0) = 
(which from equation (3.5) is independent of spin, cr); or equivalently, from equation 
(3.8): 

Et(u; = 0) = ^C/|//|. (3.9) 
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The practical importance of this condition wiU become apparent in section 5, for while 
necessary but not a priori sufficient for a SC state, its imposition as a self-consistency 
condition underlies our analysis of the SC phase, in direct parallel to our previous work 
on the normal r = Anderson model [30] . 



4. Mean field 



We start from the simplest non-trivial mean-field approximation, viz UHF as considered 
in Anderson's original paper [1]. This has two essential characteristics: that the notion 
of an impurity local moment (/x) is introduced explicitly from the outset; and that it is 
determined self-consistently, via ^ = {fii^ — ni^)^ (with (...)o an average over the mean- 
field ground state). There are three reasons for first considering this superficially simple 
one-body approximation. First, the UHF Green functions form the bare propagators 
for the dynamical many-body approach developed in sections 5ff. Second, a physical 
understanding of its limitations, often alluded to but rarely exposed, underpins what 
is required to go successfully beyond it. Finally, we show that even UHF by itself has 
virtues, producing a number of non-trivial predictions that are in qualitative accord with 
sophisticated approaches; notably that only LM states arise for r > |, in agreement with 
detailed NRG calculations [21,28]. 

The essence of UHF is that the self-energies So-('^) (equation (3.8)) are purely static: 
only the Fock term is retained, and S° = — |(Tf/|/i|. The impurity Green function at 
UHF level, Goioj) = ReGo{uj) — i7rsgn(a;)Z^o('-^), is thus 

Go(^) = i[^TH + ^iH] (4-1) 

where 



cr. 



-1 



g^{u;) =^uj + i77Sgn(a;) + -C/|/x| - A(a;)J (4.2) 
with corresponding spectral densities D^{u) = — 7r~^sgn(ci;)lm^o-(cc;) given by 

D» = [v + Mu^)U-' (43) 

[uj + ^U\ii\ - A^{uj)f + [v + Muj)f 



1 

Quite generally there are three energy scales in the problem, viz Aq"'' , U and D; we 

1 

choose to rescale in terms of Aq""", defining for later purposes a reduced interaction 
strength and bandwidth respectively by 

u = ^ (4.4, 

Al-r A 1-T- 

^0 

1 

as well as a dimensionless frequency, a) = co'/Aq"''. 



4.1. Mean-field phase boundary 

The local moment \n\ is of course found self-consistently, as described below. If = 
thereby arises, then from equations (4.2,3) the UHF propagators and spectra reduce to 
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those of the non-interacting hmit summarized in section 2.1: D^{oj) = dQ{u!). If ^ 
by contrast, equation (4.3) (with equations (2.5,9) for Aj/r) show the low-frequency 
behaviour of -D°(a;) to be 

D'^iu) (4.5) 

independent of spin a. Equation (4.5) thus gives the low-cu behaviour of the full UHF 
spectrum Dq^u) = | [i5|(c(;) + D'j'(a;)] ; and we note immediately that D{u) ~ is in 
fact the spectral hallmark of the LM regime obtained from finite-f/ NRG calculations 
[28]. We also add in passing that for finite bandwidth D, and regardless of U, there 
are always discrete (pole) contributions to the -D°(ct;) outside the band, |a;| > D. These 
are included in all specific calculations (section 6 ff), but are of little importance to the 
problem and are not discussed explicitly in what follows. 

At pure UHF level the local moment is determined self-consistently from 

= dcu [D°(w) - Dj(w)] . (4.6) 

J —oo 

Noting from equation (4.3) that the D^(u;) depend on U and solely via the 
combination 

x = ^U\f,\ (4.7) 

equation (4.6) is thus of form 

9t 

(1/^1=) ^ = (4.8) 

The UHF phase boundary Uc = U^{t) is now readily ascertained from the a; — > 
behaviour of f{x), noting that x — > may correspond either to (i) — > at some 
finite critical Uc as one might naively anticipate; or (ii) Uc — and |//| either vanishing 
or remaining finite as C/ — > 0. With 

fix) ""^^ x"" (4.9a) 

and hence from equation (4.8) 

^"-^'"-i (4.9b) 

the above possibilities are distinguished by different values of the exponent m. If m = 1 
then Uc is finite and given by 

If by contrast < m < 1, then equation (4.9b) shows Uc — and 

U^ : < m < 1 (4.10b) 

vanishes as U — > 0. Finally, if m = 0, then Uc — but 

finite : m = (4.10c) 
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Figure 1. Mean- field pliase boundary {AqU^ /U)c verus r (for wide-band limit). For 
r > ^, solely LM states occur for all U > 0. 



tends to a finite limit as C/ — > 0+. 

All three posibilities are realized in practice, as now summarized (details are given 
in the Appendix). For r < |,m = 1 results: there is a finite critical U^{r) separating 
LM states {U > U^{r), > 0) from SC states {U < f/°(r), \fi\ = 0). For r G (|, 1) by 
contrast, < m = (1— r)/r < 1; while m = for r > 1. It follows directly that for r < | 
both SC and LM states may arise, while for r > | exclusively LM states occur for all 
U > 0. This is as found in finite- C/ NRG calculations [21,28]. It is moreover specific to 
the particle-hole symmetric case under consideration: for the asymmetric case, analysis 
of the UHF equations yields a finite U^{r) even for r > |; this is again in qualitative 
agreement with NRG results [21]. 

The resultant mean-field phase boundary is shown explicitly in Figure 1 for the 

1 

wide-band limit D = oo (which depends solely on the ratio U = [//Aq"""); for later 
comparison to NRG results, the figure shows the critical AqU^/U{= f/^~^), versus r. 
For r < ^, the critical U^{r) is finite and given from equation (4.10a) using equations 
(4.8,6) and (4.3) by 

4 r J, A.M|.-AkM| ^ ^ ^^ ^^^ 



U^r) 7rJ_^ ([cu- Ar(cj)]' + A2(cc;))' 

Since Aj/r ~ |c<j|'' as cj — > 0, the low-cj behaviour of the integrand is ~ l^l"^*^; the integral 
thus converges for r < ^, and U^ij-) ~ (1 — 2r) as r ^ |— , producing the square root 
divergence in the phase boundary evident in Figure 1: AqU^ /U ~ (1 — 2r)^^^'^. 

For r > I by contrast, even the simple mean-field analysis predicts solely a LM 
phase for any U > 0. Recall however that for the non-interacting hmit U — 0, the ground 
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state is LM for r > 1 but SC for r < 1 [29] (see section 2.1). Hence for i < r < 1, there 
is a critical line Uc{r) — 'separating' SC and LM states. This underlies why low-order 
perturbation theory in U about the non-interacting limit is inapplicable for r G [|, 1] 
[29] (see also section 7.1); and the subtlety of f/ ^ is already evident at mean-field 
level by the fact that the UHF vanishes in a non-analytic fashion, viz 

1^1 ^^V^ :re(i,l). (4.12) 

For r > 1 by comparison, the ground state is LM for all > 0; and remains finite 
as f/ — > (equation (4.10c)). In fact, as shown in the Appendix, the U — 0-|- limit of 

Ia*I is 

\^ji\^q : [/ = 0+, r > 1 (4.13) 

where q (equation (2.10b)) is precisely the weight of the uj — Q pole in the non-interacting 
single-particle spectrum dQioj). 

Finally, we clarify the familiar remark that UHF is a static mean-field 
approximation. In one sense this is trivial: the self-energies T,„{uj) (equation (3.8)) are uj- 
independcnt at UHF level, being given by S° = —}^aU\ii\. But the corresponding single 
self-energy S(a;), defined conventionally via G{u}) = [uj -f- i?7sgn(a;) — A(a;) — E(a;)]~^, is 
given from equation (3.7) by 

EHF(^) = ^(^)(lf/|^|)2 (4.14) 

with g{ui) the non-interacting Green function, equation (2.3). Thus, even at UHF level, 
S(cj) is uj- dependent. For r > 1 in particular, the leading low-u behaviour of g{uj) is 
given from equations (2.3,10) by g{uj) ~ q / {uj + irisgn{uj)) . Hence, from equation (4.13), 
the leading low-a; behaviour of ^■'^■^(ti;) as U — > is given by 

E"^^ (u) ~ — : : r > 1. (4.15) 

And simple though it is, this result is not trivial: it recovers exactly the leading low-u; 
behaviour of S(a;) obtained [29] from second order perturbation theory in U about the 
non- interacting limit (which is applicable for r > 1 [29]). 



Jf..2. Deficiences 

Its virtues notwithstanding, the limitations of UHF by itself are of course severe. If 
the self-consistent mean- field local moment = 0, UHF reduces trivially to the non- 
interacting limit; there is thus no hint of the low-energy Kondo scale symptomatic of 
the r > SC phase and evident in the generalized Kondo resonance appearing in the 
modified spectral function A{u;) — \u;\'^D{u;) [29]. But the acute deficiences of simple 
mean-field arc already evident in the UHF phase boundary of Figure 1. From finite-C/ 
NRG studies [21,28] it is known that the critical AqU^ /U vanishes linearly in r as r — >^ 
i.e. Ao,c ~ r or ?7c ~ 1/r. For r = this recovers the well known fact that the normal 
Anderson model is a Fermi liquid for any non- zero hybridization strength Aq (or finite 
U), with a LM phase confined exclusively to the atomic limit, Aq = 0. This is not 
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(a) 'elastic' 



AAAAA/-> 



(b) 'spin-flip' 



Figure 2. Schematic of scattering processes as discussed in text. 



however captured by UHF, which instead produces a critical Aq/U = I/tt for r = 
(see Figure 1) and thus a spurious transition between the SC (or Fermi hquid) and LM 
phases at a finite couphng strength. 

More generally, the r — > behaviour of the NRG phase boundary for the Anderson 
model [21,28] is indicative of that for the corresponding Kondo model as considered 
originally by Withoff and Pradkin [22]; for since Aq/J/c oc r as r — > 0, the Anderson 
model can here be mapped onto a Kondo model [21,28] via the usual Schrieffer-Wolf 
transformation [32], with an exchange coupling J oc V-^i^/U. Thus, as found originally 
via poor man's scaling for the Kondo model itself [22], there exists an infrared unstable 
fixed point at Jc oc r such that for J > Jc {J < Jc) the ground state is SC(LM). Clearly, 
however, no vestige of this Kondo physics is captured at UHF level. 

The manifest deficiences of the simple mean-field approximation naturally stem 
from its static character, and in particular from the complete omission of dynamical 
spin-fiip scattering processes illustrated schematically in figure 2. Upon addition of, 
say, a J,-spin electron to an |-spin occupied impurity, two subsequent processes may 
occur, (a) The added |-spin may hop off the site, leaving behind the original f-spin. 
This is essentially elastic scattering; it is well captured by UHF alone, (b) However the 
f-spin electron originally present may also hop off the site, leaving behind a spin-flip on 
the impurity, the energy cost for the spin-fiip being on the order of the Kondo scale. This 
process, involving correlated electron motion and dynamical coupling of single-particle 
excitations to low-energy spin fiuctuations, is entirely absent at UHF level. Inclusion 
of such, to which we now turn, is however essential to circumvent all the limitations of 
the static mean-field approach outlined above, and in particular to recover the correct 
physics of the Kondo (or spin-fiuctuation) regime. On the other hand one should not 
abandon UHF entirely, but rather use it as a starting point for a dynamical many-body 
approach; use of it in this fashion is, as we shall show, necessary to ensure a successful 
description of the problem for all r > and from weak to strong coupling interaction 
strengths, U. 
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5. Dynamical self-energies 

The interaction self-energies ^^(a;) consist, as in equation (3.8), of a static Fock 
contribution (alone retained at mean-field level) plus the dynamical contribution ^^(a;) 
on which we now focus. The most important class of diagrams contributing to the Eo-(a;), 
and that we retain in practice, is shown in Figure 3a; mean-field impurity propagators 
(given by equation (4.2)) are denoted by solid lines, and the impurity interaction U by 
a wavy line. The physical content of Figure 3a is clear: having, say, added a cr-spin 
electron to a -cr-spin occupied impurity, the latter hops off the impurity, generating in 
consequence an on-site spin-flip, before returning again at a later time; and where all 
ladder interactions of the resultant particle-hole pair — reflecting the created spin-flip 
— are included. This class of diagrams thus captures the dynamical spin-flip scattering 
mentioned above (Figure 2b), and known e.g. from poor man's scaling [33] to be essential 
in describing the Kondo limit of the normal r = Anderson model. This is further 
evident from the equivalent recasting of Ti^{uj) shown in Figure 3b, which translates to 

Et(cu) = U^ g^{uj - ^i)n-+(cui) (5.1a) 



— oo 



S^(a;) = U^ r g^(u; - a;i)n+-(a;i). (5.1b) 

J-oo ^TTl 

These embody dynamical coupling of single-particle excitations to low-energy spin 
fluctuations, since the transverse spin polarization propagators U'^" /U'"^ — given as in 
Figure 3c by the ladder sum of repeated particle-hole interactions in the transverse spin 
channel — contain as will be shown the low-energy spin flip scales that are the essence 
of both the SC and LM phases. 

The transverse spin polarization propagators are given in turn by an RPA form 

with °n the bare polarization bubble (first diagram in Figure 3c). t["'"~(a;) is given 
explicitly by 

°n+-(u;) -ijy^ SiMG^iu;, - u) (5.3) 

while follows by interchanging the spin labels, t'*^i; ^'^'^ ^ simple change of 

variables in equation (5.3) gives 

°n+-(u;) = °n-+(-u;) (5.4) 

— which naturally applies also to the full n+~(a;)/n~+(a;). Only one propagator, say 
^I"'"~(a;), need thus be considered explicitly; the other follows from it. Further, since 
the real/imaginary parts of tl"' {u) are related by the Hilbert transform 



o^+_. . ^ r dc^ Im ^+-(a;i)sgn(a;i) 
y_oo TT uji-uj- i?7Sgn(a;) 



(5.5) 
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Figure 3. (a) Class of diagrams for T,rr{uj) retained in present work. Mean- field 
impurity propagators are denoted by solid lines, the on-site impurity t/ by a wavy line, 
(b) Equivalent recasting, including ingoing/outgoing propagators, to illustrate spin-flip 
scattering involved, (c) Particle- hole ladder sum in transverse spin channel; for n~+, 
spins are reversed. 



we may focus on Im °Q+ (a;). Prom equation (5.3), separating — Oai'^) + ('^) 

into retarded/advanced components and using the Hilbert transform 

c;j(.) . r a.. (6.) 

this is given by 

-Im °n+-(cu) = 9{uj) t ' dwi L'J(a;i)L>°(c;i - uj) 

Jo 

1-0 

+ e{-Lj) da;i L'J(a;i)Df(a;i -a;) > 0. (5.7) 

Finally, note that (equations (5.1)) preserve — as they must — the particle- hole 
symmetry Sj(cj) = — E|(— a;) (equations (3.5,8)); as follows using Q^{uj) = — ^j,(— a;) 
(equations (4.2), (3.4)) with equation (5.4) for the II's. In what follows we thus focus 
exclusively on E-f(a;) which, using equation (5.4), may be written as 

E^icu) ^U' U+-Mg^{u;^ + cu) (5.8) 
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and which we now consider separately for both SC and LM phases. 
5.1. SC state 

We begin with a brief overview of our approach to the SC phase, and give further details 
below and in section 6ff. As noted in section 3, equation (3.9) constitutes a necessary 
condition for a SC state to arise for [/ > 0; we show below that it is also sufficient, and 
reduces to 

T^{uj = 0) = \u\^Ji\ (5.9) 

(where Y^iu) = ReT,^{uj)). If this equation is satisfied then the generalized pinning 
condition symptomatic of a SC state, equation (2.16b), will be satisfied. And the core 
of our approach to the SC state is to enforce equation (5.9) — which refers to a single 
frequency, the Fermi level cu — — as a self consistency condition. In practice, as for the 
r = model considered hitherto [30] , this amounts to a self-consistent determination of 
the local moment, for equation (5.9) is of form 

g{U;x)^x (5.10) 

where g{U;x) = E^(a; = 0) depends explicitly on U, and upon x — ^U\ii\ (via the 
dependence of the mean-field propagators Qa{^) upon x). With a chosen approximation 
for E|(a;) — equation (5.8) in the present work — the modus operandi is clear: for given 
r and a chosen U, solve equations (5.9,10) for x and hence \fj,\; if a solution is possible 
one has a SC state, and the U above which a solution is no longer possible gives the 
critical Uc{r) for termination of the SC phase, i.e. the phase boundary between SC 
and LM states (on the assumption, indeed found in practice (section 6), that solely LM 
states arise for U > Uc{r)). 

An initial illustration of what results is seen in Figure 4 where, for r = 0.2 (and 
the wide-band limit D = oo), we show the resultant spectral density of transverse spin 

excitations, Imn+~(a;) versus u — cu/Aq^'', for U — U/A^'"" — 9, 10 and 13. The inset 
to Figure 4 shows for comparison the corresponding Im °Q+~(a;) associated with the bare 
polarization bubble, shown explicitly for U — 1?, only (since those for the other C/'s differ 
insignificantly from it). Two principal points should be noted. First, as one expects, 
Im {uo) consists simply of a high energy Stoncr band centred on ~ U\^\ ~ 10^. 
For Imll^ (cj) by contrast, it is seen that the vast majority of the spectral weight has 
been transferred to a low-cj resonance peaked at a characteristic spin-flip scale that 
is at least 3 orders of magnitude smaller in the examples shown: this is essentially the 
Kondo scale characteristic of the SC state; it will be investigated in detail in the following 
sections. Second, upon increasing U in the SC phase, uj^^ progressively diminishes, and 
vanishes at the critical C/c('^)(— 15.8) where the resonance in Imn+~(a;) becomes an 
isolated pole at a; = precisely; as discussed below and in section 5.2, the latter is the 
characteristic signature of the doubly degenerate local moment state. 
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Figure 4. Imn+-(a;) versus Cb = w/Aq'-'' for r = 0.2 and U = U/A^''' = 9, 10 and 13 
(right to left); there are no further spectral features outside range shown. The spin- flip 
scale (I'm (peak maximum) progressively diminishes upon increasing U, and dim — ^ 
as the SC/LM transition is approached. Inset: corresponding Im°n+~(a;) for U = 13, 
showing the high-energy Stoner band. Note the very different scales in the two figures. 



We now return to consider S|(a;), equation (5.8). In the SC phase, 11+ (tu) obeys 
the same Hilbert transform as ^~^~{uj), equation (5.5); using this, together with 

i_oo 27n u-uji±iri 

equation (5.8) reduces to 

/OO J\ 
Imn+-(a;) [6^01(00^ +00) + e(-u;i)g+(u;i + u)] . (5.12) 
-00 

Since Im^^(a;) = ^t:D^{u!)9{±u!) (see equation (5.6)), equation (5.12) yields 

E^(u;) = Ef (a;) - isgn(u;)Ef(a;) (5.13a) 



where 



Ej.(cu) = e{-uj)U^ f ' duji Imn+-(cui)Dj(cui +a;) 
Jo 

+ e{uj)U^ dc;i Imn+-(wi)Dj(c^i + w). (5.13b) 
Imn+~(a;) is given from equation (5.2) by 

Imn--(^) = hn^^lH ^5,4) 
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and is non-negative (as follows from equation (5.7)); so too is D^ioj) (equation (4.3)). 
Hence Ej(a;) > as required by analyticity; and E^(a;)/Sj(a;) are related by the Hilbert 
transform, equation (2.7) with F = 

To obtain the low-a; behaviour of E|-(u;) from equation (5.13b), we require that 
for -D°(cc') and Imn+~(a;). For the SC phase, UHe °n"'"~(a; = 0) < 1 (as discussed 
further below); hence from equation (5.14), Imn+~(a;) cx Im °Q+~(a;) as a; — > 0, with 
Im °Q+~(a;) given by equation (5.7). In what follows we consider explicitly the case 
where |//| > self-consistently, for three reasons, (i) This is naturally the case relevant 
to strong coupling (large U) behaviour for any r > 0; in particular to the SC/LM phase 
boundary (section 6), to the asymptotic behaviour of the Kondo scale (section 6), and 
to the consequent universal scaling behaviour of single particle spectra in the SC phase 
(section 8). (ii) The case where = (self-consistently) differs only in detail from 
> 0; the main conclusions reached below hold also for = 0. (iii) We consider the 
case of 1/^1 = in section 7.1, where we show in particular that as ?7 — > our description 
of the SC phase is perturbatively exact to (and including) second order in U about the 
non- interacting limit. 

Prom equation (5.7), using equation (4.5) for D^{ijj), the low-o; asymptotic 
behaviour of Im °Q+~(a;) is 

A H 2 



Im 



[LO] ~ 







5(r)|cu|^+'" (5.15) 

where B{r) = y/TTV{l + r)/ [2^+'^^V{l + r)] ; hence 

Imn+-(cj) oc |cj|^+^'' :cj^0. (5.16) 
The low-a; behavior of E|(a;) then follows from equation (5.13b), viz 

{t]{u) =) E}(u;) (X : a; ^ (5.17) 

with a prefactor independent of whether uj — > 0-|- or 0—; and from the Hilbert transform 
equation (2.7), 

Ef (cu) Ef (a; = 0) - 70; (5. 18a) 

with 7 = — {dYj^{uj) / dui) given by 
r°° duj El(c^) 

7 = / — >0. (5.18b) 

The full interaction self-energy S|(a;) = E^(a;) — isgn(a;)E|(a;) is given by E|(a;) = 
— |C/|//| -|- E|(a') (equation (3.8)), so the low-o; behaviour of E|(a;) is again given by 
equation (5.17); while if equation (5.9) is satisfied, the asymptotic behaviour of E^(a;) 
is given from equation (5.18) by 

Ef (u;) -70;. (5.19) 

Both Ef (cu) and Ej(a;) thus decay to zero as a; — > more rapidly than Ai/r ~ \ujY 
for any r < 1. These are the requisite conditions upon Eo-(<^) for a SC state to arise 
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for [/ > 0, as discussed in section 3. Hence, z/the basic self-consistency equation (5.9) 
admits a solution, this is a sufficient condition for a SC state. 

An important characteristic of the SC state follows directly from the above analysis, 
namely that as cu — > 

E|(a;) =E;(a;) = E(a;) : u; ^ 0. (5.20) 

That S-|-(a;) = Ej^(a;) to leading order follows directly from equations (5.17,19) using the 
basic symmetry equation (3.5); from this, using equations (3.2,3) for G{uj) together with 
the definition equation (2.12) of the single self-energy S(c<j), it follows that T,„{uj) = Ti{uj) 
as cj — > 0. Equations (5.17,19) thus give the low-c<j behaviour of S(a;), and encompass as 
a special case the Fermi liquid behaviour characteristic of the normal (r = 0) Anderson 
model, viz E^(a;) oc o;^. As will be shown in section 5.2, the above behaviour is in 
marked contrast to that characteristic of the LM phase. 



5.1.1. SC state: stability condition Before turning to the self-energy in the LM phase 
(section 5.2), and to provide continuity to that discussion, wc comment on the evolution 
in the SC phase of the local moment determined self-consistently from equation (5.9); 
and its associated implications for the low-energy spin-flip scale illustrated in Figure 
4. 

From the Hilbert transform equation (5.5) appropriate to n+~(a;) in the SC phase, 
it follows that 

Ren-(. = 0)=r^^!=5^ >0 (5.21) 

which is positive definite since Imn+~(a;) > 0. But n+~(a;) is given by equation 
(5.2) whence, since Im°n+-(a; = 0) = 0, Ren+-(a; = 0) = Re°n+-(a; = 
0)/ (l — f/Re^"*"" (cij = 0)). For the stability condition equation (5.21) to be satisfied, 
< [/Re°n+-(cu = 0) < 1 is thus required. And an explicit expression for Retl"'" {lj = 
0) is readily deduced from equation (5.3), using equation (5.6) together with the identity 
g^{uj) - Giiuj) = -U\i^\g^{uj)gi{uj); namely 

URe °n+-(a; = 0) = — / dcu [D^{uj) - £>f(a;)] (5.22a) 

(5.22b) 

where f{x) {x = ^U\n\) is thus defined, and has been introduced in section 4.1. For 
stability, > f{x) is thus required; and since f{x) may be shown to be a monotonically 
increasing function of x — ^U\ii\, saturating to 1 as x — > oo (as is physically obvious), 
the condition = 2x/U > f{x) thus amounts to 

> \iio\ (5.23a) 

for any given U in the SC phase; where |/io| is given by 

l/xol = fiWli^ol) (5.23b) 
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and from equation (4.8) is simply the mean field (UHF) local moment (denoted from 
now on by |/xo|)- 

Hence, for the SC state stability condition equation (5.21) to be satisfied, the 
local moment determined self-consistently via equation (5.9) must exceed the 
corresponding UHF moment |/xo|- This is correctly found in practice upon solution 
of equation (5.9) using equation (5.8) for E|(ci;) (as shown explicitly in Figure 6 below); 
and is reflected in turn in the spectral density of transverse spin excitations, Imn+~(a;), 
which (see Figure 4) is characterized by a strong resonance centred on the low-energy 
spin-flip scale > 0. But on increasing U in the SC phase towards the critical 
Uc{r) above which solution to equation (5.9) is no longer possible, the self-consistently 
determined approaches |;Uo| continuously from above and — ^ (Figures 4, 6). 
Ki U = Uc{r) precisely, |^| = |^o| and = 0. This is the transition point: here 
UcKe °Q"'"~(cj = 0) = 1, and the resonance in lmH^^(u;) becomes an isolated pole at 
a; = 0. The latter is the natural signature of the LM phase, since for a doubly degenerate 
LM state with flnite weight on the impurity there is no energy cost to flip a spin; it 
persists through the LM phase, where \iJ,\ — l/^oL consideration of which we now 

turn. 



5.2. LM state 

We begin by stating the result for n+~(a;) in the LM phase. It consists of a continuum 
contribution, denoted ^n+~(a;), and an a; = pole with poleweight Q > 0; speciflcally, 

n+-(a;) = ^+ ^n+-(cu) (5.24a) 

with 

dRe tl+-(a;) 



Q 



duj 



> 0. (5.24b) 



a;=OJ 

Im^n+~(a;) is given by equation (5.14), and the real/imaginary parts of ^n+~(a;) are 
again related by the Hilbert transform equation (5.5). 

The pole contribution arises from equation (5.2) for H"'"~(a;) because (i) 
URc °H"'"~(c<j = 0) = 1 in the LM phase (|//| — |//o|), and (ii) the low-o; behaviour 
of Re ''n+~(a;) is linear in lu, 

R,^,.-(.)"-^ + ,(^2i|pM^ (6.25) 



a;=0 

POO 

-oo 



where 7i {dRe (uj) / du!)^^^ = da; Im {u)sgn{u) / u"^ is readily shown to be 
positive. The low-o; behaviour of Im ^~^~{uj) is again given precisely by equation (5.15), 
and for any r > decays to zero as ^ more rapidly than uj. Hence, from equation 
(5.2), the a; — s> behaviour of H+^(a;) is —Q/u with Q given by equation (5.24b); the S- 
function part of which pole is obtained by an analytical continuation cu ^ cu + iri, which 
is unique since Q > and Imn+~(a;) > necessarily. Equation (5.24) thus results. 
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Using equation (5.24) in equation (5.8) for S-f(a;), together with equation (5.11), 
gives the basic form for E|(a;) in the LM phase: 

Et(u;) = QU'^giiu) +s Et(u;). (5.26) 

The first term, arising from the pole contribution to n+~(a;), controls the low-o; 
asymptotics of E-f(a;) as shown below. The second, ^E|(a;), is given by equation 
(5.8) with n+~ — > ^n+~, and by precisely the same argument used in section 5.1 

is given by equation (5.12) with Imll"'"" Im^n^^; hence, as in equation (5.13), 
^E|(a;) = ^S^(u;) - isgn(o;)^Ej(ci;) with ^^\{uj) > given by equation (5.13b). And 
since G^iuj) = Re^j"(u;) — isgn(u;)7rDj'(c<j)6'(— u;) it follows that E|(ci;) in its entirety is 
given by equation (5.13a); with S|-(a') > as required by analyticity, and Ej^/Ej. related 
by the Hilbert transform equation (2.7). S|(u;) is given explicitly by 

(t]{u;) T,\{u;) = 7rQU^Dl{u;)e{-u;) +^ Ef(u;) (5.27) 

the low-o; asymptotics of which we now consider. 

The a; — > behaviour of ^Ej.(a;) is obtained from equation (5.13b) in parallel 

to the corresponding analysis of section 5.1. Im^n"'"~(a;) is given by equation (5.14), 
with lm°n"'"~(a;) by equation (5.15) and Re'tl"'"" (cu) by equation (5.25); and since 
[/Re°Q+"(u; = 0) = 1, equation (5.14) thus gives lm^n"'"~(a;) oc as a; — > 0. 

Hence, using -Dj'(c^) ~ as — 0, equation (5.13b) yields 

^E\{u) oc \u\^'' u^O. (5.28) 

But since D°^{u;) oc \u;\'^ as a; — > (equation (4.5)), it is the first term in equation (5.27) 
that controls the low-o; behaviour of E|(u;) for any r > (and it is only for r > that 
a LM state arises, as will be shown in section 6). It can moreover be shown that the 
poleweight Q (equation (5.24b)) is given simply by Q = |/ioP; hence, using equation 
(4.5) the a; — behaviour of E|.(a;) = E'^(a;) in the LM phase is: 

mcu) ""^^ 4Ao\u\''e{-ou). (5.29) 

Ej^(a;) = Ej^(— a;) follows directly by symmetry and we note that in contrast to the SC 
phase (equation (5.20)), Ei-(a;) and E|(a;) do no^ therefore coincide to leading order as 
a; — s> 0, the physical significance of which will be discussed in section 8. 

The corresponding real part E'^(a;), and hence E|^(a;) = — |t/|//.o| + E^(ci;), follows 
by Hilbert transformation. Here we simply note that S|^(cj = 0) < |f/|;Uo| is found 
throughout the LM phase whence, in contrast to the SC phase (equation (5.9)), 
E?|^(a; = 0) < 0. This, together with equation (5.29), enables the low-a; behaviour 
of the full single-particle spectrum D{u)) = — 7r~^sgn(a;)lmG'(ti;) to be obtained; for 
from equations (3.2,3), D{uj) = ^ [-Dt(c^) + -D|(a;)] with D^{uj) given as a; —> by 



Da{u;) ~ ^. (5.30a) 

TT [ER(a; = 0) 
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And since = 0) = — E^(a; = 0), it follows directly using equation (5.29) that 

D{uj) ^kr- (5-30b) 



TT 



SHcu = 0) 



The characteristic low-w spectral signature of the LM phase found in NRG 
calculations [28], D{u}) oc |a;|'', is thus recovered; and we add that equation (5.30b) 
holds for the LM state regardless of whether r ^ 1. For the particular case of r > 1 
however, equation (5.30b) is readily shown to be asymptotically exact as C/ — > 0. Here, 
S|''(0) is dominated by the Fock contribution of —\U\^q\\ and as discussed in section 
4 (equation (4.13) and Appendix), |//o| ^ q a.s U — >■ 0, with q (equation (2.10b)) the 
weight of the a; = pole in the non-interacting single-particle spectrum doiuj). Hence 
as t/ ^ 0: 

^^")=^'^'V(^"""'^ :t/-0, r>l. (5.31) 

This is precisely the result obtained by us hitherto (equation (5.12) of [29]) using straight 
second-order perturbation theory in U about the non-interacting limit, which is itself 
applicable for r > 1 (but not for | < r < 1 [29]). Note that this result is not captured 
correctly at pure mean-field level alone which, from equation (4.5) with = q, differs 
by a factor of 3 (the presence of which reflects the fact that equation (5.29) for Sj(cj) as 
u; — > is independent of U). That equation (5.31) is correctly recovered as a limiting case 
of the present theory is thus a non-trivial consequence of the local moment approach. 

Finally, the low-c<j behaviour in the LM phase of the conventional single self-energy 
E(a;) — defined by equation (2.12) — may also be deduced from the above asymptotics. 
We consider explicitly r < 1 (the case r > 1 will be discussed in section 7.1). Consider 
the Hilbert transform equation (2.7) for F{u;) — F^{u;) — isgn(a;)Fi(a;). If the u; ^ 
behaviour of Fi{u!) is 

Fj(u;) a\u;\^ (5.32a) 

with — 1 < A < 1, the low-o; behaviour of F^{u) is readily shown to be 

Fr(cu) -sgn(cu)tan (f A)Fi(cu). (5.32b) 

With Fi{u) = nD{u) given as — by equation (5.30b), ReG{u) = -Fr(cj) thus 
follows directly; and from equation (2.12), S(a;) = u + i?7Sgn(a;) — A(a;) — G~^(a;) is in 
consequence given asymptotically by 

S^(a;) ""-^ sgn(a;)tan (fr)E^(a;) (5.33a) 

with 



tf{u; = 0) 



(u) ~ — ==— cos^ 



{^r)\iur ■■ r<l (5.33b) 

This divergent behaviour of S(a;) as a; — >■ in the LM phase — which as just seen is 
a direct consequence of D{uj) oc as a; — >■ — is in marked contrast both to So-(c^) 
(equation (5.29)); and to the behaviour of E(u;) in the SC phase (equations (5.17,19,20)) 
where, as befits a generalized Fermi-liquid state, S^(a;) vanishes at the Fermi level, a; = 0. 
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6. Statics 

We now consider the ramifications of the local moment approach (LMA) developed in 
the preceeding sections, beginning with 'statics'; dynamics, in the form of single-particle 
excitation spectra, will be investigated in sections 7 and 8. Specifically, we consider 
here: (i) the phase boundaries between SC and LM states, including their predicted 
scahng behaviour (section 6.3) and quantitative comparison with NRG results [21,28] 

(section 6.4). As mentioned in section 4, the problem is characterized generally by 

1 

two dimensionless material parameters: the reduced interaction strength, U — C//Aq"'', 

1 

and bandwidth D — D/Aq"''; or, equivalently, by U and U/D. We seek the critical 
Uc{r) versus r phase boundaries, as a function of U/D. (n) The evolution, and critical 
behaviour, of the central low-energy spin-flip (or Kondo) scale, a;m, that is symptomatic 
of the SC state and was discussed briefly in section 5.1 (see Figure 4); this is considered 
in section 6.1. (iii) The relationship (section 6.3) between the soft-gap Anderson model 
and the corresponding Kondo model in the strong coupling (large U) regime [21,28]. 
Particular attention will be given throughout to the behaviour of physical properties at 
low r, which is of evident importance in connecting to the 'normal' (r = 0) Anderson 
model; and for which the salient results, of which the most important are asymptotically 
exact, can be extracted analytically. 

To motivate our subsequent discussion, we begin with a phase diagram obtained 
from the LMA as detailed in section 5.1, via the limits of solution to equation (5.9) 
appropriate to the ?7 > SC state. Figure 5 shows the resultant critical AqU"^ /U{= 
U^~^) versus r, obtained strictly for the wide band limit D = oo (but, as shown in 
section 6.3, coincident in practice with that for U/D < 0.1); the inset shows the low 
r < 0.1 behaviour on an expanded scale, and the main figure includes for comparison 
the static mean- field (UHF) result obtained in section 4.1 (Figure 1). 

As seen from Figure 5, the small-r behaviour of the phase boundary differs radically 
from that of UHF, which yields a spurious SC/LM transition even for r — 0. The LMA 
phase boundary, by contrast, vanishes linearly in r as r — > 0, as found in NRG studies of 
both the Anderson [21,28] and Kondo [22,27] models; the precise form of Uc{r) as r — s> 
will be established in section 6.1. For the normal r = Anderson model in particular, 
the LMA correctly recovers (Ao/f/)c = 0: the LM phase is here confined entirely to 
the atomic limit Aq = 0, and for all Aq > the system is a normal Fermi liquid with 
a strong coupling Kondo scale whose asymptotics within the LMA will be obtained 
analytically in section 6.1 from the limit r ^ 0, and shown to coincide with that arising 
from scahng and the Bethe ansatz. Upon increasing r from zero the critical AqU'^/U 
departs quite rapidly from linearity, turns upward and then terminates at r = |. The 
latter behaviour is different from that arising at mean- field level (where the boundary 
curve diverges as r |— , Figure 1), but has been reported in an NRG study of the 
Kondo model [27]; we return to this issue when making detailed comparison with NRG 
results in section 6.4. 

One small but important point should be noted regarding the LMA phase 




Figure 5. LMA phase diagram, {Aolf^ /U)c = Ul~^ versus r (wide band limit); for 
r > i, solely LM states occur for all U > {). The mean- field phase boundary is shown 
for comparison (dashed line). Inset: r < 0.1 behaviour on expanded scale including 
the r — > Kondo asymptote (dotted) of 7rr/8 (see section 6.1, equation (6.10)). 



boundaries: the same Uc{r) is correctly obtained whether one approaches the boundary 
from the SC phase {U < Uc{r)) or the LM phase {U > Uc{r)). Recall the statement 
above that f/c(^) is obtained from the limit of solutions to equation (5.9) — T,^{uj = 
0) = — appropriate to the U > SC state; such that for U > Uc{r), S|^(0) = 
cannot be satisfied self-consistently. As mentioned in section 5.2, we find by contrast 
that S^(0) < throughout the LM phase. For the same Uc{f) to arise coming from the 
LM phase, U > Uc{r), thus requires S!^(0) — > 0— as C/ — > Uc{r)+ (on the assumption 
of continuity). This is indeed precisely as found in practice. 

The final qualitative feature of the phase diagram shown in Figure 5 is that for 
r > I solely LM states occur for all f/ > 0, as is indeed found from NRG studies 
[21,28] and seen already at mean-field level (Figure 1). Here we find, consistently, that 
S|^(ci; = 0) < for all finite U, and vanishes only as — > (where the interaction 
self-energy is of course zero by construction). 

6.1. Strong coupling asymptotics 

We now consider the asymptotic behaviour of the SC phase in the spin-fiuctuation regime 
of strong- coupling (meaning large U). From the phase diagram Figure 5, this formally 
entails consideration of the important limiting behaviour r — > (where Uc{r) oo), 
since charge fluctuations can strictly be neglected only as C/ — > oo; although in practice 
the results obtained below naturally hold over a flnite-r range. Our aim is to obtain the 
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critical IJc{t) and, relatedly, the behaviour of the low energy spin-flip or Kondo scale 
uj^{r) as [/ — > Uc{r)— from the SC phase. No restriction is imposed on the bandwidth, 
D, and the impurity level = —U/2 may lie within or outside the band; to encompass 
which we define A, used in the following, by 

A = min [D, f ] . (6.1) 

The key to extracting the strong coupling asymptotics naturally lies in S|(a;), given 
for the SC phase by equation (5.12); and in particular in E^(a; = 0), whose large-C/ form 
is readily deduced from two properties of the transverse spin polarization propagator 
Imn"'"~(c<;) (itself illustrated in Figure 4). First, that in strong coupling the spectral 
weight of Imll^ {u) is confined entirely to frequencies u > (as already evident in 
Figure 4), and in consequence the first term on the right hand side of equation (5.12) 
for E^(a; = 0) is dominant in determining its asymptotics. Specifically one finds 

— Imn+-(a;) ^ 1 (6.2) 



TT 

in strong coupling, which behaviour refiects physically the saturation of the local 
moment, |//| — > 1: it is straightforward to show, and physically rather obvious, that 
the moment saturates for Ai(a; = A)/|ei| ^1; i.e. 

I » AoA^ (6.3) 

Second, as illustrated in Figure 4, the strong resonance in Imn"'"~(c<j) occurs on the 
low-energy spin-flip scale, uj^, that diminishes rapidly with increasing U and vanishes 
as U ^ Uc{r)— from the SC phase. Imll^ (cj) is thus in practice non-zero only on 
the scale of Um', and on scales of this order KeQj^{uj) is a slowly varying function of 

from equation (5.12) by 



frequency. Hence, the strong coupling behaviour of S?(a; = 0) is given asymptotically 



Ef(u; = 0) ~[/2Re6?r(^m) — Imn+-(u;i) (6.4) 

The [/-dependence of u!ui{r) — and in consequence the critical Uc{r) — can now 
be determined from equation (5.9), which (as discussed in section 5.1) ensures that the 
generalized spectral pinning condition characteristic of the SC state, equation (2.16), is 
satisfied; for in strong coupling, where |//| 1, equation (5.9) reduces using equation 
(6.4) to 

WReG^iuJ = ^U. (6.5) 
This can be solved for uj^-^{r) once R,eQj^{uj) is known, and to which we now turn. 



Re^j^ (lo) is given from equation (5.6) by the one-sided Hilbert transform 
Reg7{u)^ [ duiD'^MP ( — - 

^ J-oo \UJ-UJi 



(6.6) 
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whose u; behaviour is dominated by that of D|(a;i) as uji — >■ 0, given generally by 
D^yiuji) ~ (Ao/vra;^)|a;i|'' (equation (4.5)) with x = ^U\fj,\] and from the latter alone it 
can be shown that the asymptotic behaviour of Re [^_p(a;) — ^j7(0)] is given exactly for 
r < 1 by 



Re [^r(-) - -^-^l-r + 0{\.\) (6.7) 



sm vrr 



— although this by itself does not give Re^^ {u = 0). An approximation to Re^j^ (a;) 
that correctly captures this asymptotic behaviour is 

Reg^iu) ^ r dcui \ui\'P ( -^—] (6.8a) 

where we have introduced a high-energy cut-off of order A = min[i5, -j]- We emphasize 
that the specific cut-off used is wholly inessential to the following arguments: the 
important point is that the prefactor to the integral is precisely Aq/ttx^. Evaluation of 
equation (6.8a) for a; > gives 

Regiiu) + 0(k|/A) (6.8b) 

use of which in equation (6.5), with x = |C/ in strong coupling (|/x| — > 1), yields 

U='-^(^-^<]. (6.9) 



(6.10a) 



TT \ r sin(7rr) 
The critical Uc{r), where Um vanishes, follows immediately from 
'AoA''\ _ Tir 

(and of course holds asymptotically as r — >■ where Uc{r) — >■ oo); equivalently, in the 
reduced units U = U/A^ and D = D/A^, 



[min(£>,|- 



u 

yielding in either case 



nr 



(6.10b) 



C/c(r) (6.11) 

Trr 

as r — > 0. The behaviour equation (6.10) is evident in Figure 5 (inset) appropriate to 
the wide band hmit, where the critical AqV^/U — U'^~^ determined by full numerical 
solution of equation (5.9) is compared to the small-r result nr/S. The latter is indeed 
seen to be asymptotically approached as r — > 0, which behaviour is reached in practice 
for r < 0.02 (a feature we compare to NRG calculations in section 6.4); wc also add in 
passing that the asymptotic region of linearity in r widens upon introduction of a finite 
bandwidth D, as will be shown in section 6.3. 
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For U < Uc{r) in the SC phase, the spin- flip or Kondo scale a;„ 
directly from equation (6.9), being given as r ^ by 

U 



hkewise follows 



(6.12a) 



= A 



ttU 



(6.12b) 



and vanishing as U - 
limit of r = yields 



Uc{r)— with the characteristic exponent 1/r for r > 0; while the 



0) = min [D, ^] exp 



TT 



U 



8Ao 



(6.13) 



which is the Kondo scale characteristic of the normal Anderson model and is 
exponentially small in strong coupling. The prefactor, A = min ^] , merely reflects 
the high-energy cut-off used in equation (6.8a); but the exponent of —ttU /SAq is exact, 
agreeing precisely with the result obtained from the Bethe ansatz [34] for the wide-band 
limit of the r = Anderson model and, more generally, with poor man's scaling [33] 
(see e.g. [2]). Note moreover that recovery of this exponent hinges on the asymptotic 
validity of equation (6.10) for the phase boundary as r ^ 0, which we likewise believe 
to be exact. 

Finally, while the results above are strictly valid as r — > 0, the exponent of 1/r for 

the vanishing of the Kondo scale cOmir) as C/ — > C/c(r)— , is found to hold generally within 

1 

the present LMA. Figure 6 shows the numerically determined (Dm = ujia/ Aq"'^ verus U 
for r = 0.2 (in the wide-band limit); and as C/ — > Uc— — 15.8, careful numerical analysis 

r ^ ^ -i5 

shows a;in ~ 1 — U/Uc ■ For U > Uc in the doubly degenerate LM phase, a;in = 
(as indicated in Figure 6); but is nonetheless 'present', giving rise to the uj = pole in 
Imn"*"" {uj) throughout the LM phase (section 5.2) and in consequence to the leading term 
QU'^Qj^ i^^) contributing to the LM self-energy S|(cj), equation (5.26) (whose existence 
is in turn responsible for e.g. the asymptotically exact result equation (5.31)). Figure 
6 also shows the evolution of the self-consistently determined local moment \fi\ as U is 
increased towards and through Uc{r), together with its mean-field counterpart |/io|, and 
illustrates the discussion of section 5.1.1; specifically that > |/io| in the SC phase, 
with — > |/Xo|+ as [/ — > C/c— . 



6.2. Connection to the Kondo model 

We comment briefly on the connections to the soft-gap Kondo model, which itself has 
been studied extensively by NRG [26,27], poor man's scaling [22] and large- mean- 
fleld methods [22] ; and to which the results of the previous section are clearly related: 
equation (6.3) is just the condition for applicability of the Schrieffer- Wolff transformation 
[32] mapping the Anderson to the Kondo model; and from equation (6.10a) is manifestly 
satisfled for r -C 1 as the SC/LM phase boundary is approached. 
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Figure 6. For r = 0.2, the spin-flip/Kondo scale Um = tOm/A.^''' versus U = U/Aq~'' 
(left-hand scale) in the vicinity of the SC/LM transition; and the local moment 
versus U (right scale, solid line), compared to the mean field \^o\ (dotted line). As 
U Uc-, (Dm ~ [1 - tJ/Ucf\ and Wm = in the LM phase, tj > tJc- 



The Kondo Hamiltonian i^K = -f^host + Hg-d consists of the host band contribution 
(equation (2.1)), with spectral density 

PhostM = C|wr^p-kl); (6.14) 

and the s-d interaction Hs-a = \ JY^k,k' '^k]'^k'iS~ + cl^Ck'i-S:^ + {cl^Ck'^ - Cfc^Cfe'i) Si^ 
with an exchange couphng constant, J, given from the Schricffcr- Wolff transformation 
in the particle-hole symmetric case by J = 8V'^/U (where a constant Vik = V is taken). 
From equations (2.5,6) and (6.14), is related to the hybridization parameter of the 
Anderson model by Aq = nV'^C, whence 



7- 

ttC [/ ■ 



(6.15) 



The critical value of Aq/C/ at the SC/LM phase boundary is however given within 
the LMA by equation (6.10a) as r — > 0; hence as r — > the critical Jc{r) in the Kondo 
model is 



Mr) 



(6.16) 



with A = min [D, This is precisely the result for the Kondo model obtained by 
Withoff and Fradkin [22] from poor man's scaling with D finite, and hence X = D. 
Likewise equation (6.12b), which with \ = D may be cast as uJ-ai{T) = D {1 — Jc{r)/ J)^ , 
recovers precisely the Kondo scale obtained by Whitoff and Fradkin via a large- A?" mean 
field treatment [22] (and denoted therein by Tq) ■ The LMA for the soft-gap Anderson 
model thus recovers precisely the correct asymptotics of the corresponding Kondo model 
in the limit r — > 0. 



A Local Moment Approach to magnetic impurities in gapless Fermi systems 31 

3 I 1 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 

r 

Figure 7. LMA phase boundaries {/^qU'^ /U)c versus r for (top to bottom) U/D 
= 100, 10, 1, 10~^ and 10~^; for U/D <C 1 common scaling occurs (bold line), as 
discussed in text. 

6.3. Phase boundaries: scaling 

We return now to the Anderson model for arbitrary r, to discuss the scahng 
characteristics of the phase boundaries predicted by the LMA, and illustrated in Figures 
7 and 8. Figure 7 shows the critical AoU'^/U{— Lf^'^) versus r curves for five values of 
U/D: 10~'^, 10~^, 1, 10 and 100; from which is seen that for U/D <C 1 (and in practice 
for U/D <0.1-1) the phase boundaries collapse to a common curve, namely the wide- 
band limit shown in Figure 5. By contrast, and for the same U/D values. Figure 8 
shows the resultant critical AqD^/U{= D^'/U) versus r curves; from which it is evident 
that foT U/D ^ 1 (and in practice U/D > 10) common scaling of the phase boundaries 

again arises. Thus, although the model contains two independent parameters — viz 

1 

U — U/Aq~'' and U/D — it is clear that the phase boundaries exhibit one-parameter 
scahng for the two distinct regimes U/D <^ 1 and » 1, according to whether the 
impurity level |ej| = U/2 lies respectively well within or outside the band. 

The behaviour found generally above is in fact suggested by the asymptotic low- 
r result of section 6.1, viz (AoA^/f/)c = iir/S with A = min[D, Moreover, the 
scaling illustrated in Figure 7 for [/ D has been observed in NRG calculations by 
Gonzalez-Buxton and Ingersent [21], who find excellent scaling for U/D = 0.2 and 
0.02 (see Figure 20 of [21], where for the symmetric model the quantity po^c there 
plotted is precisely {8/7i){AqX^ /U)c with A = U/2). These authors have also given an 
argument, based upon poor man's scaling for the soft-gap Anderson model [23], as to 
why [AqU'^ /U)c = Ul~^ should be universal; although we note that a simpler argument 
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Figure 8. LMA phase boundaries (AoD'^/U)c versus r, for (top to bottom) U/D = 
10~^, 10^^, 1, 10 and 100; ior U/D ^ 1 common scaling ocurs (bold line), as discussed 
in text. Inset: limiting U/D ^ 1 phase boundary on expanded scale, compared to 
r — > Kendo asymptote (dashed line) of wr/S 

explains this: as the bandwidth D ^ oo — which is possible for r < 1 (section 2.1) 

— this scale in effect drops out of the problem, which thus depends solely on the 

1 

dimensionless ratio U = U/Aq~'" . 

For U/D ^ 1 hj contrast, both U and U/D arc relevant in determining the scaling 
curve, but occur solely in the combination AqD^/U{= U'^'~^{D /UY). This too is readily 
understood, since for |ej| = U/2 ^ D the impurity- host coupling is controlled by the 
hybridization Ai(D) = Aq-D'' which, togther with U, sets the natural energy scales 
upon whose ratio the phase boundary thus depends (NRG calculations for the regime 
U/D ^ 1 have not to our knowledge hitherto been reported, and will be given in a 
subsequent publication [31]). We emphasize however that the limit U/D — > oo does not 
imply the total suppression of charge fluctuations close to the limiting SC/LM phase 
boundary of Figure 8, since AqD^' /U = $(r) remains in general finite at the transition. 
Charge fluctuations only become negligible as r ^ 0, where $(r) r and the condition 
equation (6.3) for applicability of the Schrieffer- Wolff transformation holds. For this 
reason the phase boundary of the Anderson model, even as U/D — > oo, differs from 
that for the Kondo model save for r — > 0. Nonetheless, as one expects physically for 
U/D ^ 1, the critical AqD^/U remains closer to its 'Kondo asymptote' of 7rr/8 over a 
wider r-range than arises for U/D <^ 1 (Figure 5). This is seen in Figure 8 where the 
Kondo asymptote is reached in practice for r < 0.1, and departure from it is modest for 
all r < i. 
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6.4.. Comparison to NRG 

We now compare quantitatively the LMA phase boundary with NRG results previously 
reported for the regime U/D <^ 1 by Bulla, Pruschke and Hewson (BPH) [28] and 
Gonzalez-Buxton and Ingersent (G-BI)[21]. BPH have obtained the phase boundary for 
U/D = 0.001, and G-BI for U/D = 0.02 and 0.2. These U/D's are sufficiently small 
that one is in the universal scaling regime where (see Figure 7) the critical AqU^ /U 
is independent of U/D; as noted above, this is demonstrated explicitly in Figure 20 of 
[21]. The NRG results for {AqU^' /U)c versus r are shown in Figure 9, together with the 
corresponding LMA result (for the wide band limit U/D ^ 0, as given in Figure 5); the 
inset shows the data for r < 0.3 on an expanded scale, and includes the r — Kondo 
asymptote of 7rr/8. It is evident that, except for the important regime r — > 0, the two 
sets of NRG data do not coincide, although they should: with increasing r the BPH 
results he increasingly above those of G-BI. We beheve however that the former results 
progressively overestimate the phase boundary with increasing r, support for which will 
be given in a subsequent paper [31]. 

From Figure 9 it is seen that the LMA phase boundary is in excellent agreement 
with the G-Bl results for r < 0.3, where the NRG points essentially lie on the LMA 
curve. Further, as discussed in section 6.1, the latter coincides in practice with the 
Kondo asymptote oinr/S for r < 0.02; and this concurs also with the two lowest r > 
points from BPH, r = 0.01 and 0.02 (see Figure 9 inset). 

For r > 0.3, the NRG phase boundaries {AqU'^ /U)c — Ul"^ increase more rapidly 
than their LMA counterpart. That such a deviation should occur is not in itself 
surprising, since our specific LMA is of course approximate and, although charge 
fluctuations are obviously included, it seeks primarily to capture the strong coupling 
physics of the spin-fluctuation regime that is asymptotically dominant for small r. As 
r |— , the LMA {AqU"^ /U)c tends to a constant value (of ~0.42, see also Figure 5); 
whereas BPH and G-Bl report a divergence in the NRG phase boundary — as indeed is 
found at mean-field level where {AqU''' /U)c is known analytically (section 4.1) to diverge 
as (1 — 2r)~2 as r — > |— , and which is also included in Figure 9. For reasons that are 
evident from the results assembled in Figure 9, we do not however beheve that the NRG 
data themselves warrant such a conclusion: to distinguish numerically between a weakly 
divergent phase boundary and a finite limit — the latter of which has been reported in 
NRG studies of both the symmetric soft-gap Kondo model [27] and the corresponding 
two-channel Kondo model [21] — is a delicate matter. We add however that this remark 
does not presume a definite answer to the question, for it is possible that the present 
LMA may not handle adequately the approach to r = |, which is without doubt a subtle 
limit worthy of further study. 
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Figure 9. For U jD ^ 1, comparison of LMA phase boundary (Ao[/''/J7)c versus r 
(solid line) with NRG results: G-BI [21] (circles), BPH [28] (squares); lines connecting 
NRG points are a guide to the eye only. The mean-field boundary, diverging as r ^ i — , 
is shown for comparison (dashed line). Inset: results for r < 0.3 on expanded scale, 
including the Kondo asymptote 7rr/8 (dashed line). 



7. Dynamics: single-particle spectra 

The ability to describe successfully single-particle spectra constitutes a stringent test of 
any approximate many-body theory. This is true even for the normal r = Anderson 
model, where as mentioned in section 1 current theories have had somewhat limited 
success; and which is but a limit of the naturally more subtle generic case of r > 0. 
In this section we consider illustrative impurity spectra, D(uS). arising from the present 
LMA, on all energy scales and for both SC and LM states. Low-frequency spectral 
characteristics in the SC phase, and in particular the predicted scaling thereof as the 
SC/LM phase boundary is approached, will be investigated in section 8. 

For a representative r < |, we begin with an overview of spectral evolution upon 



decreasing IJ — U / ^^'"^ through the LM phase, into the SC state and towards the 
weak coupling limit; including comparison of spectra on either side of the SC/LM phase 
boundary, and a brief discussion of the many-body broadening characteristic of the high- 
energy Hubbard satellites, which is correctly captured by the present theory. The weak 



limit is known to be applicable [29], we show that the present theory is perturbatively 
exact to (and including) second-order in U . Finally, comparison is made (section 7.2) 






Figure 10. Dimensionless single-particle spectra D'{Cj) = Aq D{uj) versus uj = 

oj/A^ for r = 0.2, with decreasing U in the LM phase: U = 100 (a), 25 (b), 20 (c) 
and 17 (d); the critical Uc{r) — 15.8. In (a) the corresponding mean-field spectrum is 
also shown (dashed line). Full discussion in text. 



to published NRG results for single-particle spectra in both SC and LM phases [28]. 

The procedure for determining the impurity Green function G{uj), and hence the 
single-particle spectrum D{u!) — — 7r~-'^sgn(6<;)ImG'(a;), is simply summarized: G{u!) is 
given by equations (3.2,3), with interaction self-energies f^a{^) (to which the symmetry 
equation (3.5) applies) given by equation (3.8), the dynamical contribution to which, 
E(a;), is given in the present LMA by equation (5.8); and in the SC phase the 
pinning condition equation (5.9) is enforced as described in section 5.1. In Figure 

10, for r = 0.2 (and the wide-band limit), we show the resultant dimensionless spectra 

1 1 

D'{uj) = Aq^'' D{uj) versus u = cj/Aq"'' upon progressively decreasing U in the LM 
phase: U =100 (a), 25 (b), 20 (c) and 17(d); the critical U^ir) ~ 15.8. Figure 11 
continues into the SC phase, with U =14 (a), 7 (b) and 2 (c). 

For U — 100 (Figure 10a), D'{lu) is entirely dominated by the Hubbard satellites 
centred on a; = ±C//2, with no spectral structure apparent in the vicinity of the Fermi 
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Figure 11. As for Figure 10 but for C7 < i/c in the SC phase: [7 = 14 (a), 7 (b) and 
2 (c). In (a) the corresponding mean- field spectrum is also shown (dashed line); and 
in (c) the spectrum arising from second order perturbation theroy in U (section 7.1) 
is also shown (dashed line). Full discussion in text. 



level, a; = 0. Throughout the LM phase, D'{u;) ~ \u;\^ as a; — > (see equation (5.30b)) 
and thus vanishes at the Fermi level, as evident in Figure 10. Upon decreasing U in the 
LM phase, however, a narrow low energy structure develops in the vicinity of the Fermi 
level, and becomes increasingly pronounced (Figure lOb-d) as U is decreased towards 
the LM/SC phase boundary at Uc{r). This is a precursor, in the LM phase, of the 
\u}\~'^ divergence in D{uj) characteristic of the SC phase as — 0. The latter in turn 
is evident in Figure 11 for U < t/c in the SC phase where, with further decreasing U 
the Hubbard satellites progressively lose intensity, evolving smoothly to weak spectral 
shoulders and naturally vanishing entirely as t/ ^0. 

To illustrate spectral evolution as the LM/SC phase boundary is approached, D'{uj) 
(again for r = 0.2) is shown in Figure 12 for U = 16.1 (LM) and U = 15.5 (SC), i.e. for 
U /Uc{f) — 1± S with S ~ 0.02 <^ 1; the inset compares the low-frequency behaviour 
of the spectra for |a;| < 1 x 10"^. As can be seen, the LM and SC spectra are near 
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Figure 12. Single particle spectra D'{uj) versus uj close to the SC/LM phase 
boundary. For r = 0.2, with U = 16.1 (LM, solid line) and 15.5 (SC, dashed line); 
the critical I7c(r) = 15.8. Inset: low-w behaviour on a much expanded scale. As 
5 = \1 — IJ /Uc\ — > 0, LM and SC spectra coincide to arbitrarily low energies. 



coincident on all frequency scales save the lowest (Figure 12 inset); and as 5 ^ the 
LM/SC spectra coincide to arbitrarily low energies, whence the spectra evolve smoothly 
as the phase boundary is approached. 

The behaviour described above is naturally not specific to r = 0.2, and Figure 13 
shows a corresponding spectral series for r = 0.4. The principal difference to Figures 
10,11 is that, since Uc{r — 0.4) ~ 8 < Uc{r — 0.2) (see Figure 5), the Hubbard satellites 
are less well developed in the vicinity of the SC/LM phase boundary. 

Before proceeding we comment on the strong coupling behaviour of the Hubbard 
satellite bands in D{uj). In Figure 10a we superimpose the corresponding UHF spectrum 
Dq{(jj) = i[D|(a;) + D'^{u)] with D'^^u) given by equation (4.3). This contrasts strongly 
with the full LMA spectrum D{uj): the width of the LMA Hubbard satellites is 
essentially doubled, and their peak heights halved, compared to the mean- field result. 
This additional many-body spectral broadening is well known for the normal r — 
Anderson model (see e.g. [6,30]). Its physical origins reflect both processes illustrated in 
Figure 2 whereby, having added a cr-spin electron to a — u-spin occupied impurity, either 
the added cr-spin or the —cr-spin already present may hop off the site (Figure 2a and b 
respectively). The former alone ('elastic scattering') is captured at UHF level; whereas 
the latter, involving correlated electron motion is also captured with the present LMA 
and doubles the rate of electron loss from the impurity site, thus doubling the width of 
the Hubbard satellites (with concomitant halving of their peak intensity). 

The additional many-body broadening clearly occurs only if the impurity level 
hes within the band, U/2 < D (for U/2 > D the Hubbard satellites are essentially 
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Figure 13. Spectral series D'{u}) = Aq~'' D{lu) versus a) for r = 0.4 with (top to 
bottom) i7 = 3, 6 (SC), 9, 12, 15, 18, 21 and 24 (LM); Uc{r) ~ 8.0. 



unbroadened poles). And its formal origins reside in equation (5.27) for E|(u;) 
appropriate to the LM phase: as U ^ oo, the poleweight Q (and local moment |/io|) 
naturally saturate to unity, the second term in equation (5.27) is of negligable intensity 
and Sj(c(j) ~ nU^ D^^{uj)6 {—uj); using equation (4.3) for D^{uj) thus yields Sj(a;) ~ Ai(t(;) 
for frequencies ~ — ^ appropriate to the lower Hubbard band (LHB). G^{uj) follows 

from equation (3.3) with self-energy I]|(a;) + ^tI*^)' ^^"^ ^^"^ ^ ~ ~\P 

impurity Green function itself is given as C/ ^ oo by G{u}) ~ |G|(a;) (as follows from 
equation (3.2) noting that D>^{uj) is centred on a; ~ Hence, for frequencies 

u ~ -C//2 in the LHB, 

(where we have neglected AR(a;) and T^{uj) since these merely induce a small shift in 
the satellite positions that vanishes as C/ — > oo); the analogous result for the upper 
Hubbard band follows from particle-hole symmetry, G{—uj) — —G{uj). By contrast, the 
corresponding mean- field result Gq{uj) ~ \Q'^{uj) produces 

go(^) ^-"^ ( ^u} ..^ , -.uj^-^. (7.2) 

The additional spectral broadening, and consequent halving of the sateUite peak 
intensities, is evident from equation (7.1) in comparison to its mean-field counterpart 
equation (7.2); and equation (7.1) is found to provide a numerically accurate description 
of the Hubbard satellites in strong coupling, for any r in the LM phase. 

The above qualitative behaviour is not however confined exclusively to strong 
coupling, but persists in practice throughout the LM regime and (for r < |) into the SC 
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state. This is seen in Figure 11a for r = 0.2 and [/ = 14 in the SC phase: the Hubbard 
satelhtes in D{ijj) remain centred on cu ~ ±U/2 and, in comparison to the mean-field 
Do{uj) (superimposed on the Figure), the additional many-body broadening and halving 
of the satellite peak intensities remains clearly evident. Wc add further that as r ^ 0, 
where Uc{r) ~ 1/r (see equation (6.11)), the SC phase persists to increasingly large U; 
and equation (7.1) can again be shown to hold asymptotically for the SC state, starting 
from equation (5.12) for the self-energy appropriate to the SC phase and employing a 
directly analogous argument to that given in [30] for the r — Anderson model. 

7.1. Weak coupling 

We consider now the behaviour of the spectra in weak coupling If 0, focusing 
separately on the r- regimes r < |, r > 1 and | < r < 1, and recalling that for 
r < I and r > 1 conventional perturbation theory in U about the non-interacting limit 
is known to be applicable [29]. 

We begin with r < |, where as f/ — >• the ground state is SC (see e.g. Figure 
5). Upon decreasing U in the SC phase, the local moment \n\ determined self- 
consistcntly from equation (5.9) (see section 5.1) progressively diminishes and vanishes 
at a t/o = UQ{r) (< t/°(r) (see section 4.1), as found in our previous study of the 
r — Anderson model [30]). For U < Uq, = is the sole solution; and since the 
mean-field propagators ^^(a;) (equation (4.2)) depend on U and the spin a solely in the 
combination (tC/|//|/2, both they and the polarization bubbles °Q(a;) (equation (5.3)) are 
independent of both IJ and the spin indices, and are given by the non-interacting limit 
result. In consequence the interaction self-energies So-(a;) (equation (5.8)) are likewise 
(T-indcpcndcnt, and hence coincide with the single self-energy S(c<j) defined by equation 
(2.12) (as follows using equations (3.2,3,8)). From equation (5.8) for S|(cj) = S(c<j), 
the [/-dependence of the self-energy thus arises from the explicit U"^ prcfactor thereto, 
together with the ^/-dependence of (n+^(cj) =) n(cc;) =°n(a;)/(l — U^{u)) entering 
the self-energy kernel. But for U < IIq, Il{ui) may be expanded perturbatively in U; 
and its leading term, °n(a;), when used in equation (5.8) for E(a;), recovers precisely the 
result of conventional second-order perturbation theory (SOFT) in U [29] . 

Thus, while the primary emphasis of the present LMA is naturally on the strong 
coupling behaviour dominated by the low-energy spin-flip dynamics, the resultant theory 
for r < I is also perturbatively exact to/including second order in U as U 0. This 
is seen clearly in Figure 11c for r = 0.2 and U = 2, where the full LMA spectrum 
is compared to its SOFT counterpart. We also add that it is straightforward to show 
that the stability condition equation (5.21) required for the SC state is always satisfied 
therein; and in particular that Ren+~(a; = 0) evolves smoothly as U passes through 
Uo{r), and is both positive definite and finite for all U < Uc{r) in the SC phase. In 
consequence, the LMA spectrum D{u!) likewise evolves smoothly upon decreasing IJ in 
the SC phase. 

We now turn to r > 1 where the ground state is LM for all If, including the non- 
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Figure 14. Comparison between LMA (solid line) and SOPT (dashed) spectra, for 
r = 1.5 with t/' = 7.5 X 10"^ and U/D = \; ['K{Uqf /12]D' {u) versus u) is shown on a 
logarithmic scale. As C/ — > 0, the LMA recovers SOPT precisely. Full details in text. 

interacting limit f/ = (see section 2.1 and [29]). Here too the present LMA recovers 
asymptotically the results of conventional SOPT as U — 0; as seen clearly in Figure 
14 for r = 1.5 with = 7.5 x 10"^ and U/D = |, where the LMA and SOPT spectra 
are compared. The spectra are shown on a logarithmic scale to bring out clearly the 
low-o; behaviour, which as U ^ is given by the SOPT result equation (5.31), viz 

{ApD{iu) =) D'{Cj) - [12/7r(t/g)2]|cD|'-; to which end [7r{U qf/ 12] D'{Cj) is plotted. As 
seen from Figure 14, the difference between the LMA and SOPT spectra for the chosen 
IJ is barely perceptible, and becomes even less so with further decreasing U. 

In contrast to the U ^ behaviour for r < |, however, this result is rather 
remarkable. For r < ^, with a SC ground state as U — > 0, the self-energies ^^.(a;) 
coincide precisely with each other and with the single self-energy S(cc;), as shown above: 
all self-energies are equivalent. But this is not the case for r > 1. Here, even as 
U —>■ 0, the local moment |/io| remains finite (see equation (4.13)), Il'^~{uj) is not 
therefore expandable perturbatively in U, and the self-energies S|(c<;) and ^i{uj) do 
not coincide. Nonetheless, from a knowledge solely of E|(c<;) = — [/|/Xo|/2 -|- S|(c<;), the 
single self-energy S(a;) may still be obtained directly from equation (3.7) which is quite 
general; as now considered for arbitrary interaction strength, focusing on the salient 
low-frequency behaviour. The — > behaviour of E|(a;) is given for the LM phase 
by equation (5.29), and for r > 1 it is straightforward to show that the corresponding 
real part has the low-o; behaviour E|'(u;) — S|^(0) ~ —70; with 7 = ^{11) > given by 
equation (5.18b). Using this in equation (3.7) leads to the following low-a; behaviour of 
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E(a;) = T^{uj) - \sgn{u)T}-{uj) for r > 1, 



(0) 



7rg 



(0)) ^ 



l + g7(C/) 



a; 



3?^ 



Ai(u;)' 



(1 + ?7(W 



(7.3a) 



(7.3b) 



(with corrections 0(|cj|''^^; |a;|) in the latter case); where g = [1 — ((9Ar(cj)/9cj)(^=o]~^ 
given by equation (2.10b) is the poleweight in the non-interacting single-particle 
spectrum. 

Equation (7.3) is quite general. But as C/ ^ 0, ^{U) may be shown to vanish and 
E^(a; = 0) is dominated by the Fock contribution of — |C/|/xo|, with |/io| — > as C/ — > 
(see equation (4.13)). Hence, for [/ — > 0, equation (7.3a) for example reduces to 



[7rg5(cu) + 3AogVn ■ 



(7.4) 



This is precisely the result obtained from conventional SOFT for r > 1; see equations 
(5.10,11) of reference [29]. And we note that its correct recovery reflects the fact that the 

r~ 1 ^ 2 

prefactor to the tu-dependence of Y}'^{uj) in equation (7.3), viz S|'(0) ~ [|C/|//o|] 
as [/ ^ 0, is 0(f/2) as ^ 0. 

The latter remark also sheds some light on the inapplicability of conventional 
perturbation theory in U for ^ < r < 1, since within SOFT the prefactor to the 
cij-dependence of E^/^(c(;) is 0(f/^) by construction. The low-c<j behaviour of 
in the LM phase for r < 1 has been considered in section 5.2 and is given generally 



by equation (5.33), exhibiting a characteristic 

2 



\UJ\ 



^' divergence as a; ^ with a, U- 



dependent prefactor that is again 



Ef(0) 



and for ^ < r < 1 the LM phase occurs 



for all U > (Figures 5,7,8), hence U ^ 0+ can be considered. As for r > 1, 



Sf(0) 



again as ?7 — > 0; but for ^ < r < 1 the local moment |/io| 



Ef(0) 



2r 
[/2r-l 



IS 



as 



given by equation (4.12) and itself vanishes as U ^0, whence 

[/ — > 0. For I < r < 1 the prefactor to the a;-dependence of E^/^(a;) thus vanishes as 
f/ — > 0, but with an exponent e = 2r/(2r — 1) that is in general non- integral and strictly 
greater than 2. Such behaviour cannot by construction be captured by SOFT (or indeed 
by conventional perturbation theory to any finite order in U); an inability that is not 
unexpected when one recalls that the ground state of the strict non-interacting limit 
[/ = 0, is SC for all r < 1 [29]. 

We emphasize, however, that while the behaviour in the LM phase of the single 
self-energy E(a;) — obtained as a byproduct of the LMA via equation (3.7) — is both 
singular and dependent upon whether r ^ 1 (cf equations (7.3, 5.33)), the same does 
not hold for either the self-energies central to the LMA (see section 5.2) or the 

single-particle spectra; in particular, the low-cj behaviour of D{u}) for all t/ > in the 



LM phase is D{uj) 



~ \U!\ 



(equation (5.30b)) for any r > 0. 
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7.2. Comparison to NRG 

We now compare predicted LMA spectra with the NRG results of Bulla, Pruschke and 
Hewson [28] , obtained for a fixed U/D — 10~^ using a discretization parameter A = 2 
(A — > 1 recovers the continuum hmit) and with ~ 800 states retained at each NRG 
iteration. The two sets of spectra are shown in Figure 15 as a function of cj/D, on a 
logarithmic scale to show clearly the low-a; behaviour; for (i) r = 0.25 and Aq = 0.02 
(SC state, U ~ 0.18); (ii) r = 0.25 and Aq = 0.0002 (LM state, U ^ 85); (iii) r = 0.75 
and Aq = 0.02 (LM state, U = 6.25 x 10"^). The characteristic u ^ behaviour of the 
two phases is clearly seen, Di^uj) ~ \u!\^ (LM) and ~ l^-^l^^ (SC). The two LM spectra 
have pronounced Hubbard satellites, that for r = 0.75 being centred precisely u — U/2; 
while in the SC spectrum — which is a weak coupling example {U ~ 0.18) — the 
satellites are no longer a distinct feature and have been absorbed into the band. 

The agreement between the LMA and NRG results is self-evident. For the SC 
example in particular, the agreement is essentially perfect for all uj when a small, 
controlled degree of spectral broadening is used to smooth the (necessarily discrete) 
NRG data [35]; and the ^ behaviour of the NRG data is readily shown to coincide 
precisely with that of the non-interacting limit (equation (2.11)), as we have shown is 
required for the SC state (see section 2.2 and [29]). 

Finally, we note that the present LMA includes in 11^(0;) the sum of all particle- hole 
interactions in the transverse spin channel (Figure 3a) — that captures the low-energy 
spin-flip physics (section 5) — and one can of course additionally include repeated 
particle-particle and 'bubble' interactions (see figure 9 of [30]). However as in reference 
[30] for the normal r = Anderson model, these have a very minor effect on predicted 
spectra and are not therefore considered explicitly in the present work. 

8. Spectra: pinning and scaling 

We turn now to a rather subtle, and superficially hidden, aspect of the problem: scaling 
of single-particle spectra in the SC phase (r < |) as the SC/LM phase boundary is 
approached. Such behaviour is of course well known for the normal r = Anderson 
model (see e.g. [2,4,15]) where with increasing f/ the width of the Kondo resonance in 
D{uj) becomes exponentially small, reflecting the exponential diminution of the Kondo 
scale ojk (see equation (6.13)); and the Kondo resonance becomes a universal function 

of Uj/uiK- 

The r = universal scahng curve has three essential characteristics, (i) It is pinned 
at the Fermi level: 7rAo-D(a; = 0) = 1 for all U; (ii) with characteristic low-frequency 
Fermi liquid behaviour: D{ijj) — -D(O) ~ [uj/ujkY for '^/'^k ^ 1- (iii) On larger [uj/uk)- 
scales the spectrum follows a Doniach-Sunjic law [13-15] indicative of the orthogonahty 

catastrophe, whereby Dioj) ~ [|u;|/ci;x]~" with a = 1 — 2[(5o/7r]^ and 8o = 7i/2 the Fermi 
level phase shift; i.e. D{w) ~ [|cc;|/ti;x]~2 . This has been observed in both a QMC study 
of the r = symmetric spin-^ Anderson model [15] and an NRG study of D{u!) for 
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1 0"' 1 0"^ 1 0"' 

Figure 15. Comparison between LMA spectra (Figure 15a) and NRG spectra [28] 
(Figure 15b) for fixed U/D = lO'^. (i) r = 0.25, Aq = 0.02 (SC, solid line); (ii) 
r = 0.25, Ao = 0.0002 (LM, dotted); (iii) r = 0.75, Aq = 0.02 (LM, dashed). Spectra 
are shown on a logarithmic scale, as a function oi oj/D 

o; < in the asymmetric case [14]; in practical terms, Doniach-Sunjic (DS) behaviour 
in the scaled D{u!) sets in for cu/ujk ^ 1 [14,15]. 

The question arises as to the whether the above behaviour is specific to the r — 
Anderson model, a conventional Fermi liquid; or whether it is but a particular example 
of behaviour generic to the SC phase for any r < |. We show it to be the latter, and 
to be directly and generally apparent not in D{lv) itself but in the modified spectral 
function ^(a;) = |c<;|'"Z^(c<;). That this is so is natural, since (a) it is y4(a;) that is pinned 
at the Fermi level a; = for any r > where a SC state obtains (section 2.2 and [29]); 
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Figure 16. T{lo) verus H) = cj/Aq"'' in the SC phase; for r = 0.2 (wide-band limit) 
with U = 5, 10 and 15. Note the spectral pinning (J^(0) = 1) and diminishing width 
of the Kondo resonance with increasing U towards Uc{r) — 15.8 where 0. 



and (b) in A{uj) the unrenormahzed \io\~'^ divergence in D{uj) has been removed, thus 
'exposing' directly the low-frequency many-body renormalization characteristic of the 
Kondo effect, and hence the Kondo scale as shown below. 
Specifically, we focus on T{uj) defined by 

T{uj) = ttAo [1 + tan2(|r)] \ujYD{uj) (8.1) 

which, in the SC phase, is pinned at the Fermi level (equation 2.16b): !F{io = 0) = 1 for 

1 

any r. Figure 16 illustrates versus Co — (^/Aq"'' for r = 0.2 in the wide-band hmit, 
with U = 5, 10 and 15 (< Uc{r) ~ 15.8). Spectral pinning is evident. More significantly, 
so too is a generalized Kondo resonance in J-{u) which has (a) cusp behaviour as a; — > 0, 
and (b) a width (Kondo scale) that narrows progressively upon increasing U towards 
Ucir) and vanishes as f/ ^ Ucij-) — . We show below that the latter is the spin- flip scale 
'^ra^f) discussed in detail in sections 5 and 6. 

First we estabhsh the low-o; behaviour of T{uj), and its cusp characteristics. This 
follows from D{ijj) — |[D|(a;) -|- -D|(a;)] with D„{uj) = — 7r~-^sgn(c<;)ImGo.(c<;), using 
equation (3.3) for Ga{oj) together with equations (2.5,9) for Ai/R(a;) and S^(a;) ~ —70; 
as a; — > (equation (5.19)). (Since ^^^(a;) decays to zero more rapidly than T^{uj), see 
equation (5.17), it is irrelevant to the leading low-o; behaviour.) The result is 

r{uj) 1 -sin(7rr)(9-^ + 7)|a)|^-'' - {q'^ + -ff ^{r)\Cj\^^^-'^ (8.2) 

with given by equation (2.10b) (here for r < |); and 0(r) = cos^(|r)[l — 4sin^(|r)] 
such that 0(0) = 1. For r > 0, J^{uj) as o) — > is thus dominated by an |a;p~^ cusp 
(seen in Figures 16 and 19b); but its prefactor vanishes as r — > 0, and for r = the 
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parabolic Fermi liquid behaviour J-{uj) — 1 ~ o)^ characteristic of the normal Anderson 
model is recovered. 

We now determine the behaviour of equation (8.2), as f/ — Uc{r) — 

where the low-energy spin-flip scale ujj^{r) vanishes as described in sections 5.1 and 
6.2; and for which 7 = — {d'Ef((jj)/duj) is thus required. The fuU self-energy 

S^(cj) = — |f/|/i| + T,^{uj) is given generally in the SC phase by T,^{u) = T,^{u) — S|^(0) 
(from equation (5.9)); and the low-u behaviour of T,^{uj) may be extracted analytically 
in strong coupling t/ 3> 1 where, from a trivial extension of the argument given in 
section 6.1 leading to equation (6.4), E^(a;) ~ U'^ReQ^ {uJ + uJm)- Using equation (6.8b) 
(with X = ^U) for the a; > behaviour of Re^^(t<;) then yields 



X uj^o 4Ao J. 
i:f[u;) ~ — — — rcoL 



sm 7rr 



UJ 

1 + — 



(8.3) 



(with corrections O [(a;m/A)(l + cu/cu,^)]) where A = min [D, -j] ; and hence for cu/cv,^ <^ 
1 the requisite low-o; behaviour 

tfiu) ^ -7^. (8.4) 

' sm(7rr) 

Prom equations (8.2,4) the asymptotic behaviour of J^{u;) as the SC/LM phase 
boundary is approached. If — > C/c(r) — , is thus given by 

I ~ I N l-r / . \ 2 / I ~ I \ 2(1- 



T{u) ~ 1 -4r - -— - 0(r) (8.5) 

(where the terms involving drop out since io^if) — > 0). The specific r-dependent 
coefficients here are naturally valid asymptotically as r — > 0, since the above analysis 
hold strictly as f/c(^) ^ 00 (sec section 6.1 and equation (6.11)); although in practice 
Ujyj) 3> 1 for all r < I (see e.g. Figure 5). The important point of course is that T{lo) 
exhibits universal scaling; and that, as expected physically, it is indeed the spin-flip or 
Kondo energy c<jjn(r) that sets the scale for such behaviour. Note relatedly that the 
physical content of equation (5.20) for the SC phase is the retoration, for sufficiently 
long times, of the locally broken symmetry inherent to the zeroth order mean field level 
of description; and that this timescale is l/o^m (as for the normal Anderson model [30], 
the result for which is the r = limit of equation (8.4)). For the doubly degenerate LM 
state by contrast there is naturally no such symmetry restoration, reflected in the fact 
that S|(a;) and S|(cj) do not coincide as a; — > (see equation (5.29)). 

The analysis above, while demonstrating universality, is confined to the low 
frequency behaviour uJ /u)^ <^ 1. We now consider the scaling behaviour of J-{uj) over 
the entire ur/uJm range. 



8.1. Scaling of T{u}) 

We begin with the r = Anderson model, for which !F{uj) = 7rAo-D(u;). The resultant 
J-{uj) versus 00 /uj^ obtained from from the LMA is shown in Figure 17, which universal 
form is reached in practice for t/ > 57r [30]. This scaling spectrum was itself obtained in 



A Local Moment Approach to magnetic impurities in gapless Fermi systems 



46 



1.0 



<l 



0.5 



0.0 



-30 



-20 



-10 



10 



20 



30 



co/co„ 



Figure 17. Scaling for normal (r = 0) AIM: !F{lo) = ttAoD{lo) versus Lo/iVm- The 
Doniach-Sunjic tails are seen from the (w/wm)"^ fit to the wings of the spectrum 
(dashed line). 



[30] , but one important facet of it was not noted (and wc are grateful to A E Ruckenstein 
for drawing it to our attention): the DS tails for juj^ > 1 are captured by the theory. 
This is seen in Figure 17 where an {lo juj-^'^l'^ fit is made to the wings of the spectrum. 
Detailed comparision between the LMA scahng spectrum and NRG calculations will 
also be given in a subsequent paper [31]; suffice it here to say that the agreement is 
rather good. To our knowledge the LMA is the only theoretical approach that captures 
simultaneously the characteristic low-cj Fermi liquid behaviour DiuS) — D(0) ~ [cj/ciJjJ^, 
and the DS tails for ijjjuJra ^ 1 that are known to arise experimentally (see e.g. [16]) 
but to be washed out in e.g. slave boson or 1/A^-expansion approaches [10-12]. 

We also point to the apparent small spectral feature occuring in Figure 17 at 
uj/uj-ai ~ 1. As discussed in [30] this is entirely an artifact of using the specific form 
equation (5.2) for n+~(a;) in equation (5.8) for the self-energy. This is not however an 
integral element of the LMA and (as discussed in [30]) may be circumvented, thereby 
eliminating the spectral anomaly but otherwise producing no significant effect on either 
the scaling spectrum (see figure 12 of [30]) or previously deduced asymptotics. The 
same feature is naturally present also in the r > spectra below, and can likewise be 
removed; but it is a minor effect that we are content to live with in the following. 

For r = 0.2, the resultant J^{uj) for the wide-band limit is illustrated in Figure 18 

for U — 10, 13 and 15. The inset shows the central portion of the Kondo resonance 

— whose half- width is proportional to uJui{t)-i as for the r = model [30] — versus 
1 

uj — cu/Aq"'' on an 'absolute' scale, to illustrate its rapid narrowing with increasing 
U and vanishing as C/ — > Uc{r)— ~ 15.8. The main Figure by contrast shows 
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Figure 18. Scaling behaviour of Kondo resonance for r > SC phase: for 

r = 0.2, T{uj) versus w/wm for U — IQ (dotted), 13 (dashed) and 15 (sohd). Inset: 

1 

corresponding T{uj) versus w = w/Aq"'' on an 'absolute' scale, to show narrowing of 
Kondo resonance upon increasing {/ towards SC/LM transition at Ucir) — 15.8 where 
Wm 0. 



versus uj/uj-ca-i from which universahty is evident; and note that although Uc{r) is finite 
for all r > 0, the Hubbard satellites are again eliminated from the scahng spectrum 
since 0Jm{r) — > as C/ — > Uc{r)—. We add moreover that while Uc{r) itself depends 
on the host bandwidth D (see Figures 7,8), the latter has no detectable infiuence on 
the scaling spectrum; as is expected physically, and indeed seen from equations (8.2,5) 
where the D-dependence of the former (contained in q~^) is eliminated in equation (8.5) 
as u!-[a{r) — > 0. 

Finally, Figure 19a compares the universal scaling spectra ^{uj) versus uj/u^a for 
r = 0, 0.2 and 0.4; while Figure 19b shows the data on a reduced scale u/ujj^ < 1 for five 
different r-values, to illustrate in particular the evolution of the low-o; cusp behaviour 
(equation (8.5)). From Figure 19a it is seen that, as for the r — case, DS tails again 
arise in the scaling spectrum for Lojuj^ > 1. Numerical analysis shows these to have the 
form J-{u)) ~ with exponent v{r) — \ — thus 'fiattening out' with increasing 

r as evident in Figure 19a. Note moreover from equation (8.1) that while LDiio) itself is 
not a universal function of oj /oj^, t^m-^(^) (|^|/^m)~''^(t^) does exhibit scaling. Its 
DS tail behaviour is thus ~ {\uj\/ uj^t^'^/"^ as for the normal r = model; and we note 
that this conforms to the DS law exponent of a = 1 — 2[5o/7r]^, since for r > the phase 
shift 5q precisely at the Fermi level is readily shown to be 7r/2, as for the r = case. 

For the r > SC phase generally, the local moment approach thus predicts universal 
scahng of the modified spectral function T{uj) (equation (8.1)) as the SC/LM phase 
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Figure 19. Universal scaling spectra ^{co) versus w/wm- (a) For r = 0.2 and 0.4. 
DS tails arise in all cases, and 'flatten out' with increasing r; details in text, (b) For 
cj/win < 1 and r = 0, 0.05, 0.1, 0.2 and 0.4 (from outside to inside), to illustrate 
evolution of low-w 'cusp' behaviour. 



boundary is approached; with characteristic Fermi level pinning, low-o; cusp behaviour 
and DS tails in J^{u!). In a subsequent paper [31] we shall see that these predictions are 
borne out by NRG calculations, to which detailed comparison will be made. 

9. Summary 

We have developed in this paper a many-body local moment approach to the symmetric 
soft-gap Anderson impurity model, including the 'normal' (r = 0) Anderson model as 
a particular limit [30]. The LMA is naturally non-perturbative, and both the notion of 
local moments and the a priori possibility of either a SC or LM state are introduced 
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explicitly and self- consistently from the outset; as reflected in the employment of an 
underlying two-self-energy description, together with self-consistent imposition of the 
[/-independent spectral pinning condition at the Fermi level that is characteristic of the 
SC phase [29]. 

The primary emphasis of the LMA is on single-particle dynamics — posing well 
known and hitherto unsurmounted difficulties for traditional theories — but an integral 
element of the approach also permits direct analysis of the SC/LM transition and 
associated phase boundaries. The theory offers a rather comprehensive description of 
both SC and LM phases, for any r > 0. The entire range of interaction strengths is also 
covered, including the Kondo/sp in- fluctuation physics that dominates the SC phase at 
large interaction strengths, and the 'cost-free' spin-flip physics of the LM state; as well 
as the weak coupling (small-?/) behaviour that is not as prosaic as naive expectation 
might suggest, being intrinsically non-perturbative for | < r < 1. 

While the theory leads to very good agreement with extant NRG calculations 
[21,28], a significant number of further predicitions arise from it that can likewise be 
tested, in regard both to phase boundaries and dynamics; and including for example 
the predicted universal scahng of SC spectra as the SC— >LM transition is approached 
— leading to an r-dependent family of universal spectra, of which that well known 
to arise for the normal Anderson model should represent but a particular example, 
symptomatic of generic behaviour characteristic of the SC (or 'generalized Fermi liquid') 
phase. These issues will be taken up in a subsequent publication [31], where predictions 
arising from the LMA will be shown to be remarkably well supported by benchmark 
NRG calculations. 
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Appendix 



We outline the asymptotic behaviour of the self-consistency equation (4.6) for the mean- 
field local moment as x = \U\ii\ — > 0; the SC/LM phase boundary at mean-field level 
can thereby be found, see section 4.1. The general case of a finite hybridization/host 
bandwidth D is considered, whence (section 4.1) the D^^{uj) contain pole contributions 
from outside the band (|a;| > D) for all x > 0; the corresponding poleweights are 
denoted by with +/— ioi uj > D and ui < —D respectively. 

Prom equation (4.6) the UHF self-consistency equation is thus 

= I du [Df (cu) - Dl{u)] 



D 



where (from equation (4.3)) 



J-D 



{cu + x- An{u;)y + A^.iu) {to - x - Ar(cc.))2 + Ajicu)\ 



(A. la) 

(A.lb) 
(A.lc) 



(A.2) 



It is straightforward to show that the pole contributions Qi 
we thus focus on /b(a;). 



0{x) asx ^ 0; 



(a)0<r<l 
Prom (A.2), 



dhix) 



(A.3) 



dx J, 7Tj_i,^"([u;-An{u;)]^ + Aj{u)f' 

The low-o; behaviour of the integrand in (A.3) is ~ |<^|~^'^, whence the integral converges 
for r < |. In consequence fh{x), and hence f{x), is 0(a;) as x ^ 0; i.e. the exponent in 
equation (4.9) is m = 1. Por r < | the mean-field Uc = U^{r) is thus finite and given 
by equation (4.10a). 



(b) ^<r<l 

The integral in (A.2) is controlled by its low-c<j behaviour, and we replace Ar(c<j) therein 
by its low-u form Ar(cj) ~ —sgn^u) P (r) Ao\uj\^ (equation (2.9)), where /5(r) = tan(-|r). 
Transforming the integration variable in (A.2) from a; to ^ = (Ao/a;)^/''|a;| then yields 
for r < 1 

"1 1 



/b(a;) X r 



dz 



nA^ ^0 



_(1 - p{r)z^y + (1 + (3{r)z^y + z^\ 



(A.4) 



This integral converges only for r > | whence, for | < r < 1, /b(a:) and hence j{x) have 
the a; ^ behaviour j[x) ~ x^ with exponent m = (1 — r)/r < 1. In consequence, \\x\ 
vanishes only as t/ — >• 0, see equation (4.12). 
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(c) r>l 

Prom the normalization condition upon D^{(jj), 

/O i-D 
duj D^^{uj) = 1 - / du; - [qJ"^ + Q^^\ (A.5) 

-D Jo 

Hence from (A.la,c) the full f{x) is given by 

nD 

fix) = 1 - 2 / du D^^iu) - [g^+) + Q^-\ (A.6) 
Jo 



where particle-hole symmetry is used (/^^(u;) = D^i—oo). But 0, /q dcu 

reduces to duj c?o(c<j) (with d^iijo) the band contribution to the non- interacting 
spectrum, see section 2.1); and for a; = 0, Q^^^ = Q^^ = Q{i is independent of a. 
Hence 

nD 

/(x = 0) = 1 - 2 / d^^iu) - 2Qo. (A.7) 
Jo 

But from normalization of the non-interacting spectrum doiuj) (see equation (2.10a)), 

+ 2 dcu d^icu) + 2Qo (A.8) 
Jo 

where the weight, q, of the u = pole in doiu) is given by equation (2.10b). Hence, 
for r > 1, fix = 0) = q; the exponent m in equation (4.9a) is thus m — 0, and in 
consequence the local moment 

1^1 ""-^ q (A.9) 

as follows from equation (4.8). 
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